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Abstract 

'■^ ■ This article provides an in-depth look at hadron high energy scattering by using 

gravity dual descriptions of strongly coupled gauge theories. Just like deeply inelastic 
scattering (DIS) and deeply virtual Compton scattering (DVCS) serve as clean exper- 
Q> ■ imental probes into non-perturbative internal structure of hadrons, elastic scattering 

0^ I amplitude of a hadron and a (virtual) "photon" in gravity dual can be exploited as 

I a theoretical probe. Since the scattering amplitude at sufficiently high energy (small 

Bjorken x) is dominated by parton contributions (= Pomeron contributions) even in 
, strong coupling regime, there is a chance to learn a lesson for generalized parton dis- 

tribution (GPD) by using gravity dual models. We begin with refining derivation of 
Brower-Polchinski-Strassler-Tan (BPST) Pomeron kernel in gravity dual, paying par- 
ticular attention to the role played by complex spin variable j. The BPST Pomeron 
^ I on warped spacetime consists of a Kaluza-Klein tower of 4D Pomerons with non-linear 

trajectories, and we clarify the relation between Pomeron couplings and Pomeron form 
factor. We emphasize that the saddle point value j* of the scattering amplitude in 
the complex j-plane representation is a very important concept in understanding qual- 
itative behavior of the scattering amplitude. The total Pomeron contribution to the 
scattering is decomposed into the saddle point contribution and at most a finite number 
of pole contributions, and when the pole contributions are absent (which we call saddle 
point phase), kinematical variable {q,x,t) dependence of ln(l/g) evolution and ln(l/x) 
evolution parameters 7efr. and \es. in DIS and t-slope parameter B of DVCS in HERA 
experiment are all reproduced qualitatively in gravity dual. All of these observations 
shed a new light on modeling of GPD. Straightforward application of those results to 
other hadron high energy scattering is also discussed. 



1 Introduction 



Despite plenty of data of hadron scattering, from which various quahtative features have been 
extracted, it remains difficult to derive and understand those features from the first principle, 
QCD, formulated as a perturbation theory. Gauge/gravity duality, however, can be exploited 
to study non-perturbative aspects of "hadron" in strongly coupled gauge theories. Many 
papers along this line focus on static properties of hadrons, such as mass spectra and three 
point couplings, but nothing prevents us from using gravitational dual descriptions to study 
"hadron" scattering of strongly coupled gauge theories at arbitrary energy scale [H [21 |3]. 
Hadron-hadron scattering P, total cross sections of deep inelastic scattering (DIS) [2], form 
factors of various conserved currents [H El Ej , and saturation/unitarity [3 El El EDI HI E21 E31 
[El [151 ES] are examples of non-perturbative observables that can be studied in gravitational 
dual, but potential power of gauge/gravity duality in hadron scattering is far from being fully 
exploited so far. 

Although perturbative QCD can describe evolution of parton distribution functions 
(PDF) and generalized parton distributions (GPD), initial data of the evolution cannot be 
determined by perturbation theory. Such non-perturbative initial data for PDF can be ob- 
tained from DIS experiments (and have also been studied in gravitational dual models [2]; 
see [T7] for a list of articles on DIS in gravitational dual), but GPD cannot be determined 
even from experimental data, without some theoretical modeling of non-perturbative physics. 
GPD describes parton distribution in the transverse directions [HI EHl [201 [21] ^"^^ two parton 
correlation in a hadron in general [22l [231 [211 [2S1 [211 [2Z], and hence it is an interesting object 
on its own. In this article, we take this non-skewed GPD and more generally deeply virtual 
Compton scattering (DVCS) [221 [231 [211 [23 [221 [21] amplitude at small x as examples of 
hadron scattering, and see that gravitational dual descriptions can determine how those non- 
perturbative scattering amplitudes depend on kinematical variables such as center-of-mass 
energy, momentum transfer, impact parameter and photon virtuality. Gauge/gravity dual 
also tells us how to think about various theoretical ideas that various models of GPD have 
been based on. 

High energy behavior of elastic scattering amplitude A{s,t) of two hadrons is characterized 
by poles and their residues of its partial wave amplitude A{j,t) on the complex angular 
momentum j-plane (e.g., [23 [221 [30] ) ; the poles and residues depend on momentum transfer 
t. The poles in the j-plane have often been assumed to depend linearly in t, which is 
supported by the spectrum of mesons and hadron scattering cross sections at least for finite 
range of momentum transfer t. Given the fact that the earliest version of string theory was 
born out of efforts to describe hadron scattering, it is not surprising that some successful 
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aspects of classical Regge theory are preserved in gravity dual, string theory on a warped 
background. Notable aspects of string theory on a warped spacetime, however, include i) a 
single "Regge trajectory" of string theory on 10 dimensions gives rise to a Kaluza-Klein tower 
of infinite "Regge trajectories" on 4 dimensions, and ii) those trajectories do not remain linear 
for arbitrary negative t [3]. The non-linear trajectories immediately result in non-Gaussian 
profile of GPD in the transverse directions (see section 14.21 of this article as well as |17]). 
although Gaussian profile of GPD is often assumed in phenomenological analysis. We will 
also describe how the residues of the Regge poles are determined by holographic set-up, and 
also explain how the Kaluza-Klein tower of Regge trajectories organizes itself to become a 
single contribution with a form factor in momentum transfer t < 0. 

An extra energy scale q — photon virtuality — is available in photon-hadron scattering, in 
addition to the center of mass energy W and confinement scale A. This extra parameter 
makes theoretical understanding of the non-perturbative amplitude interesting. The scatter- 
ing amplitude is dominated by a contribution from a saddle point in the complex j-plane, not 
from a pole, for sufficiently large q ^ A. The saddle point value j* depends on kinematical 
variables such as W, q and t. We find, by following this dependence of j*, instead of naively 
taking small x limit or large q^ limit, that observables characterizing scattering amplitude 
such as ln(l/g)-evolution parameter 7efr., ln(l/x)-evolution parameter Aefr. and t-slope param- 
eter B show qualitatively the same behavior in the strong coupling regime (gravity dual) as 
expected in perturbative QCD or observed in HERA DVCS experiment. As the saddle point 
value j* and the leading poles are both given by the kinematical variables of the scattering, 
crossover from the saddle-point phase to the leading pole phase may also be expected, when 
the photon virtuality decreases to a smaller valuell 

This article is organized as follows. Section |2] explains the set-up of gravitational dual for 
our calculation of photon-hadron scattering amplitude, while summarizing conventions and 
providing brief mini-reviews. Section [3] explains how the 2-body-to-2-body scattering am- 
plitudes are given in the gravitational dual set-up, and presents an explicit form of Pomeron 
kernel; this section is largely a repetition of the contents of [3], but we believe that some 
small improvements are also made and subtleties clarified in derivation and final expres- 
sion of the amplitudes and kernel. The photon-hadron scattering amplitude is discussed for 
zero skewedness in sections l4.1H4.3j section 14.11 explains momentum-transfer t dependence 
of the imaginary part of the scattering amplitude, while its impact parameter b dependence 
(i.e., transverse profile) is described in section Hl2l Section [4.31 is devoted to the real part 



^ A similar crossover behavior has already been observed in the real part to imaginary part ratio of the 
hadron- virtual "photon" scattering amplitude in gravity dual We will elaborate more also on this 

crossover behavior in section 14.31 
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of the amplitude. The interpretation of the scattering amphtude on a warped spacetime in 
Froissart-Gribov-Regge language on 4 dimensions is given in detail in section I4.1.3t while 
section 14.1.41 is devoted to the t-slope parameter of the scattering amplitude. We will see in 
section 113] that GPD can be identified within the scattering amplitude even in the strong cou- 
pling regime, and discuss the form factor of GPD. In section [5l we address a question whether 
there is anything we can learn about GPD of the real world from the GPD calculation in 
the strong coupling regime. Section |H] briefly describes straightforward application of various 
results and observations in this article to other high-energy hadron scattering processes. 

2 Model 

Such gravitational backgrounds as [211 [321 [331 [31] are an ideal framework for holographic 
calculation of certain types of hadron interactions, as they have built-in confinement mech- 
anism at IR, and allow renormalization group interpretation in terms of holographic radius. 
Those models have background geometry that are approximately AdS^ x W with some com- 
pact 5-dimensional manifold W, apart from the region near the IR boundary. The hard wall 
model |2] replaces these geometries with an AdS^ x W backgrouncj^ 

ds^ = GuNdx^'dx^ = Qmndx'^dx'' + {gw) abdO'' dO^ (1) 

d2 

gmndx'^dx'' = e^^'^'\r,^^dx>'dx'' + dz^), e^^^'^ = — . (2) 

An IR boundary is introduced at ^ = A^^ and boundary conditions on various fields are 
set by hand instead. The dilaton vev is simply assumed to be constant, e*^ = Qs- Such a 
background is not obtained as a stable solution to the Type IIB string theory, but Type IIB 
string calculations on such a background (while ignoring NS-NS tadpoles) are expected to 
maintain qualitative aspects of certain hadronic processes in the original holographic models. 
We use the hard wall model in the rest of this article, as it makes it possible to compute 
various physical quantities and study dynamics without consuming too much time. 

Deeply virtual Compton scattering (DVCS) and double deeply virtual Compton scattering 
(DDVCS) are elastic scattering of a hadron h and a (virtual) photon, 7*(q'i) + h{pi) — > 
'l^*H'l2) + h(p2), Figure [H with the kinematics 

ql > A^, g2 = (DVCS), ql ql > A^ (DDVCS). (3) 

As the target hadron, we use one of normalizable modes in a scalar degree of freedom 0(x, z) 

^ In our convention, 77^^ = diag(— 1, 1, 1, 1). We use M,N,... in labeling coordinates of ten dimen- 
sional spacetime, m,n, . . . for the AdS^ part, and fj^v, p,<7, k . . . for the coordinates of the 3+1 dimensional 
Minkowski spacetime. 6°-, 9'', - ■ ■ are dimensionless coordinates of W, and {gw) abdd"" dO^ is the metric of W . 
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Figure 1: a cartoon picture showing elastic scattering of DVCS / DDVCS, with momentum 
labels on the external lines 



on AdS^ in the holographic set up. In particular, we use a scalar field z) originating 
from ten dimensional SUGRA fields upon Kaluza-Klein reduction on W, just like in [35]. In 
case of reduction of dilaton, for example, 



(l){x,z,e) = (l){x,z)Y{9) 



(4) 



with a non-trivial harmonic function Y{6) on W, the target hadron h corresponds to a 
glueball|§ 

The holographic wavefunction of the target hadron h is obtained by solving equation of 
motion derived from an effective action of the scalar field (f){x, z) on AdS^. The bilinear part 
is 



Srh. 



d'^xdzy/^ 



2^1 



where 



R^Yo\{W) 



O 



and a dimensionless constant^ cs of order unity and a mass parameter M'^ depend on the 



K 



IIB 



Nl 

i?3 



(5) 



(6) 



^ To study DVCS and DDVCS of baryons, Z?-branes should be used in the holographic set up, instead of 
a Kaluza-Klein mode of Type IIB SUCRA fields. We hope that there is still something to learn in an easier 
study with a scalar glueball target. 

^ It is defined as follows: 



(7) 
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choice of Y{6). When a Dirichlet condition^ is imposed on (f){x^z) at the IR boundary 
z = A~^, the wavefunctioij§ of a hadron hn (corresponding to the n-th normahzable mode) 
with incoming momentum is given by 

^{x, z) -> e--$.(.), V5^$„(.) = ^^"' ^Vlo-ritr 



where J^{x) is the Bessel function, the n-th zero point of J^{x), and A = 2 + a/4 + M2i?2. 
Mass of the hadron /i„ is given by 

= AjA-2,n- (9) 

We will not specify the excitation level n of the target hadron hn{x), as we will pay attention 
only to qualitative aspects of DVCS / DDVCS amplitudes, not to numerical details that 
depend on the excitation level n. 

As the (virtual) photon probe of DVCS / DDVCS amplitudes, we gauge an _R-symmetry 
associated with W, and use it as a probe, just like in [2]. The type IIB metric field on 10 
dimensions become a massless vector field on AdS^ through 

SGma{x, Z, 9) = Am{x, z)Va{0), (10) 

where v"'{6)d/d6"- is a Killing vector of W, and Va = R^{g{w))abV^ ■ The kinetic term of the 
effective action of Am{x, z) on AdS^ is given by 



R^cbCa 



Sa = — — / d'^xdzy/^ 



2Ki 



,lp prnn 



(11) 

with dimensionless coefficient^ ca and cb- The non-normalizable wavefunction of the vector 



^ Qualitative aspects of our results will not change, when a Neumann boundary condition is imposed, 
instead. 

^The wavefunction $„(z) satisfies dz-y/— 5e^^"^$„(z)$,„(z) = (5„m, so that the hadron field hn{x) 

in (/)(a;, z) = /i„(a;)$„(z) has a canonical kinetic term upon dimensional reduction to 3+1 dimensions. 

^ Arbitrary chosen normalization of the Killing vectors Vi = vf{d/d9"-) does not remain in ca, as we define 
it by 

T^hrnit^t,) CA = ——^-== [ d^'OV^^vy'^gwiO); (12) 

the generators t^'s are chosen so that they satisfy the same commutation relation [ti, tj] = ifijtk as the Killing 
vectors, {vi, Vj} = fijVk, and the Cartan metric remain the same. In the case a non-Abelian subalgebra of the 
Killing vectors of W is gauged, FmnF"^^ in the kinetic term above should be understood as T^^trfl;(Fmn^™") 
using the generators we explained above. 



5 



field Am{x, z) for the "photon probe" with incoming spacehke momentum is given by 

^Fp,{x,z) ^ zcj(g,e^-e,g^)gz|KiM + ^^|;^/i(gz)|e^'^-^ (13) 

^ cj{qh, - {q ■ e)q,)z I^Ko{qz) - fj^jj^^ ^oiqz)^ e^"-, (14) 

where cj is a dimensionless constant of order unity, q = \/q^, and e(g)^ is the polarization 
within the four dimensions. For sufficiently spacelike much larger than A^, the second 
terms Ii{qz) and Io{qz) are negligible in (fT3l [T4l) . but the full expression needs to be used for 
the final state on-shell photon in the DVCS. 

The elastic scattering amplitude A('y*h — )■ 7^*''/i) of a hadron h and a (virtual) photon 
probe 7^*) is calculated in holographic model by using the world-sheet non-linear sigma model 
with the background metric([Tl [2]) and inserting vertex operators whose Born-Oppenheimer 
approximation [36], [37] are specified by the wavefunctions (jH [8]) and (fTOl [T3l [Hj) [3]. The 
Compton tensocl T^,, is defined by removing polarization vectors, 

Ai^*h-,^^*'>h) = e^Me^,r, (16) 

and the goal of this article is to determine five independent structure functionJ§ Vi 2,3,4,5 (e.g., 
[42]) in 

T^-" =V,P[q^rP[q2]; + V2{p- P[qi]np ■ Pfe])^ + ^3(^2 ■ P[gi])^(gi • Pfe])'^ 

+ V,{p ■ P[qi]nqi ■ P[q2]r + V,{q2 ■ P[gi])^(p ■ P^ + Ae^''^'' q^.q^.. (17) 

The last term with the coefficient function A should vanish for a scalar target hadron /i in a 
parity-invariant theory. Here, we introduced a convenient notation 

(18) 

Those structure functions, Vi, 2,3,4,5(2:, 77, t, g^), should be expressed in terms of Lorentz invari- 
ant kinematical variables x,ri,t and g^, where 

P" = ^(K+P^), 9" = ^K + g2"), A^ = pi^-pt = q>i-q^, (19) 

®In the real world QCD with QED probe, where only fermion partons are charged under the probe, 
{2nfSipi+q,-p2-q2)iT'''' = -J d*X2 J d^xie-'"''^' e'"'-''' {h{p2)\T{r ix2).r ixi)}\hip,)) . (15) 

^Leading order perturbative QCD result in terms of non-perturbative GPD is found in [35] in the case of 
a scalar target hadron, and in 09] in the case of a fermion target hadron. See also [JO] [41] . 



P[q] 
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x=-^, r/ = -^, s = W' = -{p + qf, t = -A^ (20) 

2p ■ q 2p ■ q 

just as in standard literature in perturbative QCD. The parameter t is assumed to be small 

(21) 

throughout this article. The two conditions on the kinematical variables ([3]) and ( |2T1) com- 
bined is sometimes referred to as generalized Bjorken regime. In this article, we focus on 
non-skewed DDVCS {rj = 0). More general case (77 7^ 0) including DVCS (77 = x) re- 
quires farther analysis. Holographic calculation of the DDVCS amplitude should reproduce 
the pure forward amplitude (whose imaginary part is the deeply inelastic scattering (DIS) 
amplitude studied in 0) when the skewedness rj and momentum transfer t are set zero; 
Im Vi{x,r],t,q'^) -> Fi{x,q^) and (g^/(2x)) x Im V2(x, 77, t, g^) F2{x,q^). 

For simplicity, we will study the sphere amplitude contribution to the four closed string 
external states. Although the sphere amplitude alone cannot discuss how the unitarity is 
maintained in the scattering, sphere amplitude is sufficient for large enough N^. (or for not 
too large s = W"^, for sufficiently large g^, or for sufficiently large impact parameter b). It is 
also possible to discuss with the sphere amplitude how the scattering approaches unitarity 
limit [11]. The pion cloud [l3lllllll5] contribution to the impact parameter dependent profile 
[16| W7\ , however, can be studied only in a more realistic holographic model containing pion 
[Ml SHI EDI EI] by examining x = 2 — 2g — h = —1 amplitude. 

3 Pomeron Contribution to the Structure Functions 

Holographic QCD has already been used to study DDVCS / DVCS amplitudes in the litera- 
ture. Reference [12] calculated DDVCS / DVCS amplitudes for scalar target hadron (a set-up 
that we explained in the previous section), with a single supergravity field in the s-channel 
resonance (see Figure [2]). This contribution to the amplitude corresponds to the resonant 
contribution to the deeply inelastic scattering in the t = limit, and is known to be domi- 
nant in the large limit for moderately large 't Hooft coupling A and moderate x [2]. We 
can certainly learn the structure functions of DDCVS / DVCS processes in this way for the 
kinematical range mentioned above, but this sugra resonance contribution in the holographic 
calculation corresponds to the higher twist double trace operators of weakly coupled gauge 
theories, and does not tell us much about the non-perturbative input of GPDs. 

The operators that are approximately twist-2 in perturbative gauge theories, on the other 
hand, have large corrections to the anomalous dimensions in strongly coupled regime (except 
the spin 2 operator). This "twist-2" contribution still dominates in the DIS amplitude for 
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Figure 2: Contribution to the DDVCS amplitude from double trace protected operators {s- 
channel sugra field resonances). (cj2x/g^)~^Im V2 is shown in (a) as functions of a; G [t], 1] 
for four different values of skewedness 77. The solid line is for t] = 0, and dashed lines are for 
?7 = 0.1, 0.2 and 0.3. We set t = in panel (a). On the other hand, the panel (b) shows the 
(— t) dependence, while we set 77 = 0. The solid line is for {—t)/q^ = 0, and three dashed 
lines are for {—t)/q^ = 0.1,0.2 and 0.3, and are shown for the range < x < (1 +t/(4g^))~^. 
g^/A^ = 10^ is used for both (a) and (b). As a target hadron, the lowest normalizable mode 
of a scalar with A = 5 was used. 



given at sufficiently small x [2], and are naturally expected to be so in the off- forward 
DDVCS amplitude as well. In this article, we study this "twist-2" contribution in the small 
X region by using holographic calculation, to get some hint on the non-perturbative input of 
the GPDs. 

The "twist-2" operators of gauge theories correspond to string states in graviton trajec- 
tory in holographic descriptions [521 13]. Contributions that involve such string states in the 
graviton trajectory to hadron scattering processes are called Pomeron amplitudes, and are 
known to be expressed as in fl62l) [2], [3l [11], [12]. In this section, we begin with refining the 
derivation of ( 162|) . Our derivation largely follows the ones in sections 2 and 3 of [3] and 
combine various improvements already made in [TTl [12], but we believe that the following 
presentation also made a couple of small improvements. All of the following issues are closely 
related: 

a) choice of integration contour in complex angular momentum plane, 

b) validity of deformation of the integration contour, 

c) origin of signature factor [1 + e~'"^^]/ sin(7rj), 

d) absence of non-sense poles at negative integer angular momenta j, and 

e) the fact that the sphere amplitude of string theory is at best interpreted as a sum of 
t-channel and u-channel exchange of particles with various spins, not purely a sum of 
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t-channel amplitudes. 

It will be made clear how we should thinl<F°l about these issues a)-e) in writing down the 
Pomeron amplitude (l62i) : interpretation provided in section 4.3 of this article is affected by 
the understanding on the issue e) we obtain here. 

Another improvement is to replace A2, a derivative operator on AdS^ often used in 
Pomeron propagator in the literature, by Aj(t) in ( 147|) for complex angular momentum 
j. The Pomeron wavefunction "^llJ for spin j G C is also introduced. Although this change 



does not leave a practical impact on the expression to be used as the Pomeron kernel (|6l 
this extra conceptual clarification of the role played by complex angular momentum j in the 
Pomeron kernel will enable us to provide clear theoretical understanding of the form factor 
associated with Pomeron-hadron-hadron coupling in section 14.11 



In section 13.21 we focus on DDVCS amplitude and find explicit expressions of Pomeron 
contributions to the five independent structure functions. It may be an option for busy and 
practical readers with some familiarity to [3l [TTl [12] to skip section 13.11 for the first reading. 

3.1 Pomeron Kernel from Sphere Amplitude 

Reference [3] derives Pomeron kernel in its section 2 by modifying Virasoro-Schapiro scatter- 
ing amplitude of closed string on 10-dimensional fiat spacetime so that it is understood as an 
amplitude of scattering on a curved spacetime background (with small curvature a' /R^ <^ 1); 
we will combine it with the discussion in section 3 of [3] to fill a small gap in the process of 
modifying the scattering amplitude on flat 10-dimensions to that on the curved spacetime. 

In fiat 10-dimensional spacetime, the sphere level scattering amplitude of two NS-NS 
closed strings has a factorized form 

A(sio,tio)(27r)i05i°(pi + gi -P2 - ^2) = K G, (22) 

where sio and tio are Mandelstam variables in 10-dimensions. The factor G is a function of 
Sio and tio, and independent of polarizations of the external strings; 

Gi^io.ho)--^ 11 r(l + a'e/4)- 

The factor on the other hand, is given by wavefunctions — momentum and polarization — of 
states involved in the scattering, whose explicit form is found in textbooks [S3]0 Normal- 
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By no means, we consider that the logical derivation in the following is the only possible one. 
We follow [2], however, to use the factor K that is s^^ times that of |53], and further includes delta 
function of momentum conservation; s'^q can be factored out from the factor K [S3], because we are interested 
only in the large sio small j^iol scattering, as in [2j[3]. 
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ization of the factor K is 



1 



[a'^gsc'^^-^] [a^gse'^n K'^^.e-^^-^ [a'^^.e-^^-ls^o (24) 



for dilation-dilaton scattering, up to a constant of order unity. The factor K for the case of 
our interest, dilaton-graviton scattering as a holographic model of DDVCS, is a little more 
complicated, because (dimensionless) polarizations of graviton external states are involved. 
We will first present the derivation of the scattering amplitude fl62p and explicit expression 
of Pomeron kernel flMl) . using the case of dilaton-dilaton scattering. Polarization dependent 
statements are deferred to section [3l2l 

We will study small x ^ g^/s DDVCS amplitude, with momentum transfer (— t) ^ g^, 
(12T]) . that is not necessarily smaller than the hadronic scale A^; here, s = and t are 
Mandelstam variables of 4D kinematics of DDVCS. This means, as we explain later, that we 
need to examine the Virasoro-Schapiro amplitude in the kinematical region a'sio ^ 1, and 
(— tio)/sio ^ 1, but not necessaril}0 |a'tio| ^ 1- Ignoring terms that are suppressed by 
(a'sio) or (sioAio)) one finds that 

r(-a'^/4) TT 1 



^ n r(l + a'^/4) sin (7ra'tio/4) r2 (1 + a'tio/4) 



TTa'tio\ , ^ / 7ra'sio\ . f na'tio 
cos ■ — + cot ■ — sm 



. (25) 



X 



In this expression, cot [ ^°'^^" ) contains s-channel poles of strings. As long ag^i < arg sio < 
TT, in |a'sio| 00, 

cot [-^j (26) 

As a result, one arrives at a well-known expression of the Regge behavior of the Virasoro- 
Schapiro amplitude on flat 10- dimensional spacetime: 

4 sm(7ra'tio/2) r^(l + a'tio/4) \ 4 J 



^^This is a difference from the application to DIS in [2], where a'tio 0(1/\/A) can be ignored. 

This is consistent with the well-known ie prescription in the quantum field theory, which defines the 
physical amplitude in the limit of argsio — >■ +0. Moreover one can understand (pS)) as the average of 

cot in a'sio > 1 as in [2]. 
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(a) 



(b) 



Figure 3: Integration contours Ci, C2 and C3, and various singularities that appear in the 
complex j-plane. The contour C2 comes from Cj — zoo, avoids the j = a(tio) singularity to the 
left, and goes to Cj + ioo; Cj G M is the real part value of the contour of C2 in the asymptotic 
region (both ends). In order for the integral ^'•'"^''(sio, tio; to return a finite value, 

we need to take > —1. 



There are a number of merits in seeing amplitudes in complex spin j'-plane in classical 
Regge theory [281 1291 [30], as well as in perturbative QCD [SU |29], and furthermore, the 
j'-plane description also has a couple of extra advantages in the present context. For one, 
we can clarify subtle points in how the derivation of Pomeron amplitude based on vertex 
operator OPE in section 3 of |3] is related to the somewhat heuristic modification of the flat 
spacetime amplitude in finding a scattering amplitude on a curved spacetime in section 2 of 
[3], as we explain shortly. Furthermore, we can find a clear guiding principle in the heuristic 
modification process of the scattering amplitude, although this process has not been crystal 
clear so far (in our eyes) in the literature. 

Let us first see, with the j-plane description, that ^(sio, tio) can be decomposed into three 
pieces. First, note that the amplitude fl27p is the same as 

where a(tio) = 2 + a'tio/2, which is the trajectory of graviton, and the contour of integration 
in complex j-plane Ci is set as in Figure [3l Mathematically, one can add to the integrand a 
function of j holomorphic at a(tio). In order to see the explicit relation with vertex operator 
OPE, we choose the following expression. 
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We will explain their relation soon, along with the meaning of f2, which is a real positive 
number. Secondly, suppose for now that 



arg(sioj 



/2. 



(30) 



Then the integrand of fl29l) vanishes fast enough when |j| goes to infinity while satisfying 
Re j > — 1. It is, thus, possible to change the contour from Ci to C3 — C2; see Figure [31 
The r(j/2) factor in the denominator is crucial in the convergence of integration at large real 
positive j; the condition ( !30|) and the choice Cj > —1 of C2 are necessary for the convergence 
of integral along C2. Finally, let us split the integral along C3 into two pieces, one with the 
Qi-a(tio) factor in the integrand replaced by 1, and the other with the same factor replaced 
by [fi-'""^*!") "-'-]• As a whole, ^(sio,tio) is given by a sum of three pieces, 

= -g^'''\s,o,ho;n^-''^'''^) + G^'^'Hs.o^ho; l) + e;(^^)(5io,tio; ^^^-"^'^"^ - 1). (31) 



Physical meaning of the second term above will be clear. Replacing the j-plane integral 
along C3 by residues at the poles j = 0, 2, 4, ■ ■ • , 



1 



j=0,2,4,- 



a{t 
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2 / / \ j-2 



1 



(32) 



which is regarded as t-channel exchange of spin j string on the graviton trajectory]^ the 
l/(j — a(tio)) factor is understood as the propagator of a string with mass = 4(j — 2) /a', 
and this contribution to the total amplitude K G is proportional to (sio)-', as the dominant 
contribution to the factor K in sio S> |tio| is proportional to s^q. 

Vertex-operator OPE in [3] in world-sheet calculation also yields the same expression as 
). Let us use (0, 0)-picture vertex operators 



vr^(o,o) 



: e 



(1) 

MiNi 



a 



a 



idX^' + — (gi ■^)^^i 
2 



:,(33) 



,(2) 
'M2N2 



a' 



N2 



-iq2-X{0) . 



(34) 



for a massless incoming NS-NS string with (^1)10 = and a massless outgoing NS-NS string 



with (^2) 
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0. The OPE of these two vertex operators includes the following series of 



^''Although the factor — a{tw)) of j = term in the summation seems to give rise to a(iio) = 
tachyon pole (i.e., a'tio ~ —4), the amphtude (sio, tio! 1) actually does not have such a tachyon pole. 



This is because a factor pyri 



ajtia) 



cancels the a(tio) = pole in j = term. 
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(a) (b) 

Figure 4: A sphere amplitude describing scattering of four NS-NS string massless states. 
Later on, the vertex operators V3 and V4 are used for incoming and outgoing hadron states, 
and Vi and V2 used for incoming virtual photon and outgoing virtual/real photon for DDVCS, 
respectively. V2,3,4 are fixed at w = 0,1, 00 on the g = world sheet, and the sphere 
amplitude is obtained by integrating over the complex coordinate {w, w) of Vi over the 
sphere. Contribution from integration over the < < 1 region (a) and the one over the 
1 < \w\ < 00 region (b) can be regarded as t-channel and w-channel amplitudes, respectively. 



operators 



.(1) . .(2)^ 



a(t 



X \w 



e'-' ■ e'-' ) I — I I — I \w\ ^ + ■■ 

a'tio/2 . ^iqi-(wdX+wdX) ^i(q-i-q2)-X g-^ . 



(35) 



In the first line, we omitted terms which are less singular than \w\ ^ in the w, ii; — )■ limit 
Keeping only terms for simplicity, we obtain. 



Vr^(u;,zz;)Vr^(0,0)~(6«-e(2)) 



,(0,0) 



a{t 



lOj 



5^-^i.^-^«;^-^O,,,(0,0), (36) 
- k\k\ 

k,k>0 

a{tio) 
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The 



terms become those with the greatest power of sio, that is, s^q ^" . Less singular terms of 

~'^w~^ . . . , have lower power of siq. This leading term in the OPE 
of string scattering in 10 dimensions becomes the leading term in small x (large s = W'^) in the structure 



OPE, whose coefficients are I w| ^,\w\ "w 



function Vi in 4 dimensions in section [3.21 Leading terms of other structure functions in small x, however, 
are determined in section without using such an OPE on world sheet. 
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when 
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Lo - 1 = - j4 + - 2, Lo - 1 = - j4 + A: - 2, (37) 



are weights of operator C^^, and 

O,,^(0,0) =: (zgi ■ aX)'= [iq^ ■ dX)' e'^''^-'^'>'' {0,0) : . (38) 

The sphere amphtude is given by integrating w over the entire complex plane. Only the 
Lq = Lq terms remain, and 



" {Vf''\w) Vf^'\0) vl-''-'\l) vl-^'-'V)) X |oo|^ (39) 

x(O,/2,,/2(0)V3(l)V4(oo))x|cX)|^ 



Because (0^/2^72(0) V3(l) V4(oo)) x |oo|^ ~ (a'sio)-', we find that the \w\ < 1 contribution 
to this world-sheet amplitude has the structure of ^'•'"^^(siq, tio; 1) in (15^ . 

Now VL in ( !28l) has a clear meaning: the cut-off of integration of \w\. The third term in 
(13T|) . ^("-^^^(sio, tio; — 1), corresponds to the 1 <\w\<VL contribution of the world- 

sheet amplitude above. Obviously the cut-off VL should be taken to infinity. Thus, the sphere 
amphtude using V1-V2 OPE should be proportional to 



G{sroMo)= lim [e^^^^H^io, ^lo; 1 - 0) + e^^^^H^io, ^lo; - 1)] • (41) 

There is no contradiction between (14T!) and (I3T1) : whenever 

- 1 < «(tio), (42) 

the contour C2 can be chosen as a straight line from Cj — zoo to Cj + ioo for some — 1 < Cj < 
a(tio). Now, it is easy to see that 

lim 6;(^^)(sio,tio;l]^'^"(*^«)) = 0, (43) 



^^We adopt a convention of |53| . where Lq a-nd Lq vanish on physical states. Thus, Lq and Lq here 
corresponds io Lq — \ and io ~ 1 in [S]- 
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because of the factor0 Although each piece of (131 p is well-defined and such relations 

as ( l3T| |4T| |43|) can be established only under the conditions (!30l l42l) . the expression (l28l) is 
always well-defined, and it should be regarded as analytic continuation of the world-sheet 
amplitudJ^ pij) off the kinematical constraints fl30l H2l) . 

We have so far discussed Regge limit of closed string scattering amplitude in flat 10- 
dimensional spacetime, but our true interest is in scattering in a curved spacetime ~ AdS^ x 
W. The prescription of section 2 of [3] is to rewrite the factor K in fl2^ by using wavefunctions 
of four external states in the scattering, and replacing sio and tio in ^(sio, tio) by appropriate 
differential operators acting on the wavefunctions. We also follow this line of argument, while 
making a couple of improvements. 

As the expression for K becomes a little messy for dilaton-graviton scattering (holographic 
dual model of DDVCS amplitudes), we postpone working out the K factor for this case until 
section [3^ For an easier case, the K factor for dilaton-dilaton scattering (holographic dual 
of elastic glueball scattering) becomes 

2^113 J 

= (27r)^5^(pi + qi-p2- q2) 



dz^/^) / d'>e./^){e~'^sf{^{z)Y{e)f{^'{z)Y'{e))\ (45) 



where (^{z)Y{9) and (^'{z)Y'{9) are normalizable wavefunctions of the initial state hadrons, 
and Cs is a dimensionless constant of order unity, sio in (124|) is replaced by (e~^^s), including 
the warped metric!^ 

In order to obtain the factor ^(sio, tio) for scattering in the curved spacetime AdS^ x W , 
Sio is replaced by the Minkowski part (e~^^s) with the warped metric, just as above for the 
K factor. As for tio in ^(sio,tio), on the other hand, one cannot drop such terms as d1 or 
c?g, unless \t\ ^ A^ [3]. Although it is not immediately obvious which derivative operators 
should be used in curved spacetime, we adopt a prescription to go to j-plane description 



This also means that the u-channel contribution ^J'-'^^-' (sio, tio! — 1) has a finite —J- +oo Umit. 

Moreover, because the total amplitude is also equal to 

e^^^nsio,tio; 1) = ^^'^■^'(sio,tio; 1) - G^'^'-\sw,tw; i), (44) 

one can see that the u-channel (1 < |ti;| < co) contribution limo^oo [^''^^■' (•^lo, ^lo; ^^"'~"^*"''')] is equal to 
—Q'^'~^^\siQ,tn)] 1), at least mathematically. The t-channel and it-channel contributions are given by integral 
along C3 and — C2, respectively, and the total amplitude is given by Ci. 

-'^^The same story is in between the Virasoro-Schapiro amplitude ([23|) and / dPw \w\~"' — w\~"' ^1"^ . 

^^Although sio also contains derivatives in z and Q, such terms are ignored, because they are suppressed 
relatively by of order A^/s and are negligible in the high-energy scattering. 
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I^T]) . first, and then replace tio in the spin j partial wave by 

tio^ A,(t) + i?"2v^ + /2-2<5,, (46) 

where 



A,{t) = e-^H + e^\-gY^l^dA-g)'l^e-^^d^e-^^, (47) 



Since the l/(j — a(tio)) factor in fl32|) is regarded as a propagator, the first two terms in fH6|) 

operators appearing in the equation of motion of spin j modes on 



are the rank-2 differentia 



the curved background!^ In addition to this reasoning based on local field theory intuition, 
it is also possible to determine j — a(tio)(= 2Lo) in direct calculation in worldsheet nonlinear 
sigma model. The last term R~^^j is a mass correction; 6j is a constant of order unity, and 
does not contain a derivative. The above prescription also can be applied to f l28|l . because 
p8|) is equivalent to ()29|). 

Combining both the factors K and Q, we obtain 

Ais,t) [ dz^d'9^{e-'^sfmz)Y{9)fg{e-'^s,t,oW{z)Y'i9)f, (49) 



1 f f a'7r\ 1 + 1 



V 4 ; sinvrj 12(^/2) 



j-a(A,(t)+i?-2(5,. + V2^))- 



(50) 



Note that all the derivatives in tio are placed on the right of all the ^-dependent factors 
y/—g and (e~^'^s)^ in the first line, and {e~'^^s)^~'^ in the third line; Aj(t) simply acts on 
the wavefunctions {^'{z)Y'{9)y. This is not asymmetric in [^{z)Y{9)f' and [^'{z)Y'{9)]'^, 
because the z derivative part of Aj(t) is Hermitian under the measure dzy^—g[z)e~'^^'^''^\ 

We will drop in the following; this is because on a compact space W has non- 
positive discrete spectrum. Apart from the constant mode, whose eigenvalue of is zero, 
eigenvalues are negative and at least of order unity. Thus, at high-energy scattering. 

In (s/A^) > VX, (51) 

a factor 

(a'e-2^s) ^ e ivr (52) 



20 rz 



-1/z and rj;^ -(5!^/z under the metric 
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in [a'e '^^s)'^ is always suppressed, unless z is extremely small. The condition (I3T]) corre- 
sponds to exponentially small x in DIS/DDVCS, 



,n(^^j»v^. (53) 

which we will assume in later sections. The scattering amplitude now has an effective de- 
scription in 5 dimensions, 

where and c^/ are defined in ([7]). 

The (a'e~^"^s)-' oc (Q!'e~^^s)"'^-''*^*^/^ factor is non-local, and one can rewrite the factor Q 
as a diffusion kernel by inserting a complete system of the operator Aj (t) |3]. Eigenfunctions 
of Ai(t), 

A.W^!?(^,^) = -^*.?(^,^), (55) 

are given by 



t(v^/A) V^-.-(v^M) 

(56) 



where /^(x) is the modified Bessel function. Here, we imposed a Dirichlet boundary condition 
at IR in order to keep the expressions a little simpler, but essence will not change when 
a Neumann condition is imposed instead]^ For neg ative@ t, the eigenfunctions ^'•^^ for 
eigenvalues —[y^ + 4)/i?^ with < z/ G M form a complete system. Replacing [$'(2)]^ by 
/ dz'S{z — 2;') [$'(2')]^ and using a relation 

di^^Sit.^'^^it,^) = [{-9Y/'e-'^^]-'5iz-z'), (58) 







The reflection coefficient R{i^,t) defined by 

{t,z)<xh, iV~tz) + R{iy, t)I-,, {V~tz) (57) 

becomes R{v,t) = — [-fiiy(Co)/^w(^o)]|^g=^/3t/A under the Dirichlet boundary condition. It satisfies R{—v,t) ~ 
\/R{v,t) and {1 + R{—v,t)) = for iv e Z. (If both \/~^ and v are real valued, R{v,t)* = R{—i',t), although 
we do not use this property in this article.) When the Neumann boundary condition dz[e^^^'i'\''J' {t, z)] = 
is imposed instead [3], R{v,t) defined as above also has all of these properties. 

We understand that is real positive for negative i, that is, arg = 0. In its analytic continuation 
to real positive t through the upper half-plane in complex then, arg ^/--t = —t: 12. 
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we find that 



4 y sinyrj T^{j/2) \ 4 



J - (j. + 5,/(2VA)) 
where e^' = e^^^'-*, and 



^ 4'!;'(«,^)*E'(t,z')[t'(^')l', (59) 



^ 2 - (60) 

With a little more work, the expression above can be made completely symmetric for the two 
initial/final state hadrons; with a notation 

5 = e-^-^'s, (61) 

one finally arrives at the expression we are familiar with in [31 [TTl [T2] : 



Cs_ IT 

where 



A{s,t) ^ I dz^^)P{z) / dz'^^^)P{z')lC{s,t-z,z'), (62) 



P(z) = c^($(^))^ P\z')=c,.{^\z')f (63) 

for glueball-glueball elastic scattering. The Pomeron kernel ]C{s,t; z, z') becomes 

1 f ,. 1 + e-^^^' 1 



JC{s,t;z,z') = -8RVX [ du—f dj 

Jo 27rz Jcuu) 



CiH sinTTj r2(j72) 



^■^(^)vl/£)(t, ^) e-^^(^')<ii^it, z'). (64) 



Because e-J'^(^)*£^(t, z) = e-2^(^)*SJ(t, ^) (see ([56])), [e'^'^^^^] in the kernel can be replaced 
by [e~^'^\l/-^^] mathematically. We should remark that this expression fl64p is valid only in 
zero-skewedness scattering. 

It should be emphasized that we have [e~'''^\l/£^(t, z)][e~^"^'^l{^{t, z')] in the Pomeron ker- 
nel, not [e-^^-^li\t,z)][e-^^-^[i\t,z')]* that uses complex conjugate. This is a big difference, 
because the integrand in (!64l) can be regarded as a holomorphic function of (j, u) (except 



18 



some pole loci), and the kernel itself is expressed as a holomorphic integral in the spin- 
anomalous dimension (j, 7) plane 1^ note that 7 = A — j — 2 = ^// — j [3]. The kernel (!64l) is 
also holomorphic in the momentum transfer t of hadron scattering, and we will exploit this 
nature in later sections. 

The j = 2,4, ■■ ■ poles from 1/ sin(-7rj) correspond to the spin j particle exchange, and 
contribute to the real part of the kernel. The so called nonsense poles j = 0, —2, —4, ■ ■ ■ of 
1/ sin(7rj) in Regge theory are canceled]^ and even become zeros in the j-plane due to the 
factor l/r^(j/2) in the j-plane representation of the kernel. The absence of nonsense poles 
can easily be traced back to the Virasoro-Schapiro amplitude fl23p and its Regge limit f l27p . 
One can also trace the origin of l/r^(j/2) factor in the argument of vertex operator OPE; 
see (|36|). 

The remaining singularity in the j-plane comes from — — Sj/{2\/X)). Sj vanishes 
at j = 2 |3], because massless graviton does not receive mass correction even in curved 
spacetime (put another way, energy momentum tensor has vanishing anomalous dimension), 
jj, is a function of and hence the zero locus of the denominator determines the relation 
between j and u. The contour Ci{u) of j integration is around this (i^-dependent) pole, 
and after integration, this pole locus determines the large s (high energy) behavior of the 
scattering amplitude. 

The pole locus jr{i^) in the j-plane is given approximately by 

Jrii^)=J. + OiX-') (65) 

for \u\ < C(l), and hence for - 2| < 0{1/V\). For |z/| ~ \^/*, and hence for ~ 0{1), 
one can still find the pole locus jr{i^) recursively [T3] : 

Jril^) = + + O(A-l) = + 0(A-V2). (gg) 

This means that the pole locus jr(^) is shifted from due to the 6j correction, not by 
an amount as large as the leading 0{1) term, but by of order (9(A~^/^). The 0{l/\/X) 
corrections are just as important as the term in Finally, for ^ A^/^, and hence 

for ^ it even becomes impossible to try to find the pole locus j = jr{v) recursively. 

It is also known that for j > -\/A, 7-j relation changes drastically [52]. To recap, the high- 
energy behavior of the Pomeron kernel is described fairly well by s^" for |z/| ^ (9(1), but the 
exponent jr{v) begins to deviate from quadratic in i/, when \v\ becomes comparable to 

^•^This is quite common in perturbative QCD; e.g., [55l [56l [29l [3] . 

^"'Thc l/r^(j72) factor also renders the total t-channel exchange contribution ^J*-'^^-' (sio, iio! 1) finite, and 
makes the decomposition possible, as we have already mentioned. 
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A^/^. This is equivalent to 

\ 1 



R / r (67) 
R J a' 

loosely speaking, that is when R, "the Kaluza-Klein momentum of Pomeron in the holo- 
graphic radius direction" becomes comparable to the string scale @ 

3.2 Extracting Structure Functions of DDVCS Amplitudes 

Although the Pomeron kernel (164|) is universal for all the hadron scattering processes at 
high energy (small x), the impact factors P{z) and P'{z') should be chosen for individual 
processes. P{z) in dH |63]) can be used for elastic scattering of two scalar glueballs [3l [T2] . 
Two independent structure functions of DIS cross section are also expressed as in (imaginary 
part of) ([n2D, using ([HI (HSD for the target hadron impact factor P'{z'); the impact factors 
P{z) for virtual photon have also been determined for the two structure functions In the 
case of DDVCS amplitudes for a scalar target hadron, there are five independent structure 
functions, and we need to determine the impact factors P{z) for these individual structure 
functions. 

As a holographic model of DDVCS scattering, we use graviton-dilaton scattering, as 
we have announced in section [21 In order to determine the impact factors for the DDVCS 
scattering, it is easiest to start from the known form of the factor K in f[22]) for the graviton- 
dilaton scattering, and carry out the process corresponding to f l45p . the same strategy as in 
|2]. For DDVCS, 

~ (27r)V(pi + gi-p2-g2) 

I dz^g I rfWV^^;.(0)^;'^(^)F^(^)F2^™(^){(Pl)p(p2). + (P2)p(Pl).}($(^)>^(^))^ 

(68) 



2^IIB 



where c'^ is a constant of order unity; Fi and F2 should be understood as wavefunctions 
in (11 3[ [T^ except the plane wave part already taken into account in the four momentum 
conservation. For the Fi for incoming virtual photon and F2 for outgoing virtual photon, 
the wavefunction with = (qi)^ and the one with = (— ^2)^ should be used, respectively. 
We have kept only the terms where four momenta of graviton (gp)'s and dilaton (pp)'s are 
contracted in (168!) . because such terms dominate in high-energy scattering, as in [2]. pi and 
P2 are made symmetric in {■ ■ ■ } in (1^5]) . because the polarization of the graviton propagator 



Another condition |sio/iio| ^ Ij which was used already at ([25]), is satisfied for i/ smaller than or 
saturating ([67| : this is because |a'tio| remains |Q!'Aj(<)| ~ {v^ + 4)/-\/A < 0(1), while we consider a'siQ ^ 
[s/A2]/VA»e^/\/A» 1. 
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is symmetric. Inserting the Regge limit process-universal amplitude Q to the right of \^~g 
and {e~'^'^q ■ p)^, we find that 

e-2V''e-^V^{(Pi)p(P2). + {p2)M)AG{e-'''sMo){<^{z)Y{e)f. (69) 

For DDVCS at exponentially small x = —q^/{2q ■ p) ~ —(f/iv + = <f/s, fl53|) . only 
the constant mode on W contributes, and the expression above is reduced to an effective 
description on 4+1 dimensions, just like in (EU CgC^c^i is now replaced by d^CAC^- 

We are now ready to read off the impact factors for the DDVCS amplitude. The impact 
factor for the target hadron is the same as in the elastic scattering of the target hadrons, and 
P'{z') = Phh{z') is the same as (163|) . The impact factor for the virtual photon is 

Piz) = CA R's-'g^-{F,r^{zmUz){{p,),{p2)^ + ip2)MU- (70) 

It can be decomposed into the five structure functions of f[T71) : the result is as follows: 

2 

V, ^ -^h, V, (71) 

qz qz 

where we treated all of (A^/g^ = —t/q^) and x to be much smaller than unity@ Here, 

J,(x,r/,t,g2) ^ j dz^^^ j dz'^^^)P%.{z) ]C{s,t,z,z') P,,{z'), (72) 

with 

P^!^.(z) = c]R'e-'^^^\{q,z){K,iq,zmq2z)K,iq2z)], (73) 

P?A^) = ' ^2 Mz){Ko{qizm<ll^)M<l2z)]. (74) 

When one ignores terms that are suppressed by powers of x, as we have done so far, the 
five structure functions are in fact given by only two contributions, Iq and Ji. Although 
the longitudinal and transverse polarization of photon in the incoming 7*(q'i) + h{pi) beam 
axis is different from those in the outgoing 7*(g2) + ^(^2) beam axis in the presence of non- 
zero momentum transfer, they become approximately the same in small angle (and hence in 
small x) scattering; the two beam axes are precisely the same in the t = limit. In this 
sense, Jq and Ji correspond to the amplitude of the virtual photons with "longitudinal" and 
"transverse" polarizations, respectively. 

^^We have already assumed 77 = 0. 
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4 DDVCS Amplitude and GPD at Small x in Gravity 
Dual 



Now that concrete expressions are given to the (Pomeron contribution to the) structure 
functions of DDVCS amphtude fl7T| [72]) . with an exphcit expression for the Pomeron kernel 
( I64p . let us evaluate the integrals to get physics out of them. The momentum transfer t 
dependencJ^ of DDVCS amplitude at small momentum transfer < — t < is highly 
non-perturbative information, and this is where gauge/gravity dual can play an important 
role. 

We will first focus on scattering amplitude (i.e., structure functions) in sections l4.1H4.3j 
the imaginary part of the amplitude is studied in sections |4. 1 1 and [4. 2[ which sheds a light on 
non-perturbative form of the parton distribution in the transverse direction at small x |21] . 
The real part of the amplitude is described in section 14.31 We will argue in section 14.41 that 
GPD can be defined as an inverse Mellin transformation of operator matrix element and is 
calculable even in gravity dual; the structure functions and the scattering amplitude as a 
whole are given in the convolution form involving GPD, just like in the QCD factorization 
formula. 

It is useful in studying the DDVCS amplitude in the generalized Bjorken regime (|3l [2T]l 
to write down the Pomeron kernel explicitly as follows. Carrying out j integral around the 
pole j = jrii^), and using the explicit form of \E'£''(t,2;) in (156|) for the hard wall model, we 
obtain 



/C(s, t] z, z) = iy/Xe 



2A-2A' 



sinh nu 



tz)L 



du 



tz' 



1 + e- 



sm vrjy 



1 



t/A) 



t/A) 



tz)Ii^ 



tz' 



(75) 



AVXe 



-2A-2A' 



du 



1 + e-''^^-^ 



sm 71 ji, 



r2(jW2) V 4 



TT 



K- 

~tz') - 



-t/A) 



hu{V^t/A) 



tz') 



(76) 



Although both of the second lines of ( 1751 1761) are equivalent to (1771) when integrated over 
those in ( |75l [76l) . which is not symmetric under the exchange of z and 2', turn out to 
be a little more convenient than the z z' symmetric expression (177|1 in evaluating the 



For very large momentum transfer —t ^ in the QCD in the real world, however, perturbative QCD 
can be used to argue rough scaling behavior [57l [SH [SH |60l [61] of the DDVCS amplitude. 
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DDVCS amplitudej^i dominant contribution to the amplitudes comes from small z region 
because of the virtual photon wavefunctions localized toward UV boundary, 
and for such a small value of (a/^2;), the Iiiy{y/^z) and I-iy{^/^z) terms in ^!'fj(t^ z) have 
quite different behavior as a function of {iu). Iiu{y/—iz) decreases rapidly toward positive 
Re (iz/), while I_iy{-\f^z) toward negative Re (iz/). This is why the I^iy{\f^z) term has 
been turned into Ii^{y/^z) in ( 175| [76]) by relabeling —lyhyu. 

4.1 Momentum Transfer Dependence of the Imaginary Part 

The imaginary parts of the structure functions simply come from the imaginary parts of the 
integrals Ji and Iq fl72|) . and their imaginary parts come from the imaginary part of [l + e"*'^-'] 
in the Pomeron kernel. It must be straightforward to substitute the expression of the kernel 
above into fl72|) and evaluate them for kinematical variables of our interest; References [21 [11] 
(see the reference list of [T7] for other articles) have done that for purely forward case t = 0, 
and we can just carry out a similar procedure of calculation for < —t case as well. Before 
doing so, however, we find it worthwhile to write down the amplitude in the following form, 
which leads us to a better theoretical understanding of the t 7^ amplitude. 

Exploiting the kinematical constraint of the generalized Bjorken regime f[2Tl) . we expand 
liuiV—iz) in the Pomeron kernel (176]) in a power series and keep only the first term. Contri- 
butions to li from the higher order terms are suppressed by powers of (— t/g^), because the 
integration over the holographic coordinate z is dominated by the region z < 1/q. Ignoring 
the higher order terms is like dropping higher twist contributions in perturbative QCD. One 
can then see that the structure functions can be written as 



/oo 
-00 



Sm TTJy 
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The second line of (|75l [76]) can also be written as 



21/ 



UuW —t/ A) 



(77) 



The second term vanishes in the A — > limit for fixed z, z' , s, t, and the Pomeron kernel on AdS^ in |10j is 
reproduced; the A — >■ limit of the Pomeron kernel in this article is different from the one in [TU] only by a 
factor of (4\/A)^~Vr^(j72), which becomes 1 at j = 2. 
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where 



A 



hh 



J dz^-g{z)P^'}^,{z)e 



K 



.2 



dz'^/-g{z')Phh{z') 



-2A{z) 

■ R 



,{j.-2)A{z') 
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(79) 
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We will see in section HX^ that CfJ and correspond to OPE coefficient and matrix 
element of a "twist-2" spin j = operatoi^ 

The two factors C\]j and A^^ depend on kinematical variables and x as follows: 



■iii) 



— c, 



^(0 



where •ju = ii^ — ju, and 



A 



^ 4 j {RAr gU^t/A) ^ 4 



X 



{qTmrati^t/A)- (82) 



ci} is a constant of order unity that depends only on v (when r/ = 0, so that q2 = qi — q), 
while the t dependence of the structure functions Jj now remains only in a dimensionless 
factor g^^{^/^/A) whose definition can be read out from (IHO j [821) . Therefore, the imaginary 
part of the structure functions are given by 



Im /, 



oo J2 



P(jV2)V4v^x; \q) ^^^^^ 



t/A). 



(83) 



4.1.1 Small momentum transfer: —t < 



Now let us evaluate the amplitudes (1831) . ffist, for the case with small momentum transfer 
< — t < A^. The purely forward amplitude (i.e., one for the DIS total cross section) is 
a part of this story. The v integration can be evaluated by the saddle point method for 
exponentially small x fl53l) . just like in [11]. The factor g^^{^/^/A) has a non-zero finite 



Because of the factor 1/Kiv{y/^R) in the integrand of (|80|) . the second hne in ([80]) is normahzed to 
be unity at UV boundary z' = R i). A careful reader may notice that Re (w) > is assumed in 

the expression of (|80| . This assumption gives rise to no problem, since throughout this article we always 
estimate the u integral (|78l) in the lower half i^-plane. 
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(dimensionless and 0{1)) limit when —t 0, and g^i^{\/^/A) is a slowly changing function 
of (ii^), unless <^ —t. Thus, large u dependence comes only from 

ji^ / A \ ^'^ 

(84) 



^l^*{q/A,x,-t<A') = VX^^^f^, (86) 

In ' 



.a/Ax 

in ( l83l) . The saddle point u* is determined by the kinematical variablejf§ x and as in 

ln(g/_A) 

'q/l 

and the amplitudes approximately become 

This leading order expression of the saddle point approximation (187|) can be improved by going 
to higher order; those terms would give rise to corrections that are suppressed relatively by 
powers of yA/ln[(g/A)/(v^x)]. 

Ignoring all the factors of order unity and the Gaussian measure of the saddle point 
approximation, we find that the DDVCS amplitude for small momentum transfer is roughly 
of the form 

Im/i(a;,r/ = 0,-t < A2,g2) _ |^__J j^_J = (^^J (x) e ^mc^/avv^.). (88) 

These results can also be used for the DIS structure functions, with -Fi(x, g^) = Im /i/2, and 
^2(0;, g^) = Im [x(/o + /i)] in the purely forward limit t = and 77 = 0. 

To characterize the dependence and x dependence of the DDVCS and DIS structure 
functions, let us introduce effective exponents, as often done in phenomenological analysis of 
structure functions. 

2n <91n[x /i(x,r/ = 0,t,g2)] . , . 2^ 91n[x/i(x, r/ = 0, t, g^)] 
^'='^-^"'''^) = ' ^-.(-,t,g) = . (89) 



As we discussed around the expression of the kernel ([75)1 is vahd for \v\ < A^/"*. Therefore, the 

kinematical variables are restricted within the following region: 

|..1<AV^^A-V^ln(^)>ln(|). (85) 
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From f l88p . we find that 



7efr.(a;,t,g^) = 7,.., \es.{x,t, q'^) = jy* - 1. (90) 

Both 7eff. and Aefr., and hence the and ln(l/x) evolutions, are controlled by the saddle 
point value iu*] the saddle point value iu* G M>o in (jHSD becomes large for large and 
decreases to zero for smaller x. 

The effective anomalous dimension 7efj. = iv* — jy* ~ {iv* — 2) for a given q^ is positive 
for larger x, and negative for smaller x. This means that the (generalized) parton density 
decreases in q^ evolution for larger x and increases for smaller x. For a given value of x, 7ce. 
turns from negative to positive as q^ increases, and the (generalized) parton density at that 
value of X begins to decrease; the parton splitting from x to smaller x' becomes faster than 
the splitting from larger x' to x. This is precisely the behavior expected in [2] . Note that the 
essence of seeing this expected behavior is in keeping the (g, x) dependence of the saddle point 
value z/* in small x (ln(l/x) ^ \/~\) and large q^ (g^ ^ A^) region, without naively taking 

— 00 limit. Parton picture still remains even in the strong coupling regime, although 
the parton contributions do not dominate in the DIS structure functions at moderate x 
(a/a < x), and the DGLAP evolution is very fast. 

Similarly, for a fixed value of g^, Aefr. becomes smaller than 1 for sufficiently small x, 
rendering the x-integration for j = 2 moment converg ent at X = [n]@ For a given value of 
X, Acff. = ju* — 1 increases for larger g^, implying that the (generalized) parton distributions 
rise more steeply toward x = for higher g^ (see Figure [H] (a)). This observation has already 
been made in the purely forward t = case [T7] . 

The discussion so far clearly shows the conceptual importance of the saddle point value 
of iu* and hence of j^* . The Pomeron amplitude corresponds to a sum of contributions from 
various states/operators of spin j G 2N, and the sum can be expressed as a holomorphic 
integral in the complex z/-plane. Equivalently, this integral can also be expressed in the 
complex j-plane, via the relation j = j^. The saddle point value j*(= ju*) reflects the "center 
of weight" of contributions from various j G 2N. The whole amplitude approximately shows 
the x-evolution and g^-evolution of "spin j* operator", as clearly shown in f l9UI) . We will see 
later in this article that not just the and g^ dependence of jes. and Aefr. in (ISOj) but that 
of almost all the observable parameters (t-slope parameter in section 14.1.41 and real part to 
imaginary part ratio in section 14. 3p of the photon-hadron scattering amplitude are governed 
by the saddle point value j*. 

■^^ One can also see the convergence directly, when the scattering amplitude is given by a complex j-plane 
integral; the a;-integration for j ^ n moments converge for all n larger than the largest real part of the 
singularities in the j-plane. 
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The essence of the saddle point approximation is in the following expression 

(91) 



27rz V / 

where we already assume that we are interested in the small x and large region, and ignored 
various factors the are irrelevant to the determination of the saddle point. The anomalous 
dimension 7(j) as a function of complex spin variable j is [3] 

7(j) = 7,^. = IV, - J = {2^/\{j - jo)}^/' - J (92) 

in the hard wall model, where Vj is the inverse function of j = jv We obtain an approximation 

2m \q) ' dj .^^^ ln(g/A) ^ ' 

The saddle point value of j is a functioijf^ of [ln(g/A)]/[ln(l/a;)], just like in fl86l) and in f l93|) . 
The saddle point approximation at the leading order ignores terms that are suppressed by 
1/ ln(l/a;), but keeps full [ln(g/A)]/[ln(l/x)] dependence at all order. Higher order corrections 
in the saddle point approximation take account of l/ln(l/x) suppressed corrections. 

One will notice that this argument does not rely on detailed form of the anomalous 
dimension 7(j) very much. Indeed, exactly the same line of argument has been used in 
perturbative QCD for the study of behavior of PDF in the small x and large region 
[5^ . The anomalous dimensions 7(j) of the twist-2 series of operators in weak coupling 
gauge theories (with some variations in the approximation scheme (e.g., double leading log 
approximation (DLLA))) are not the same as those in gravity dual models such as (1921) . but 
they can be continuously deformed from one to the other by changing the value of the 't 
Hooft coupling [3]. Discussion so far makes it clear i) that contribution associated with the 
"twist-2" series of operators (parton contribution) does exist in the weak coupling and strong 
coupling regimes alike, ii) that the kinematical variable dependence of the parton contribution 
in the small x and large region can be captured by the saddle point approximation fl9T|[93l) . 
and iii) that the (g, x) evolution of the parton contribution remains qualitatively the same 
in the both regimes, despite the difference in the anomalous dimensions. This similarity of 
the parton component of a hadron in the both regimes is an encouraging factor in trying to 
take advantage of gravity dual descriptions to study non-perturbative aspects of partons in 
a hadron of a confining gauge theory (like the real world QCD). We will elaborate more on 
this in section [51 

A little more careful discussion should be given in non-conformal theories. It will not be difficult, however, 
to incorporate the running coupling effect within weakly coupled gauge theories or separately within gravity 
dual descriptions. 
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4.1.2 Large momentum transfer: ^ —t 

Let us now evaluate the holographic DDVCS amplitudes fl83l) for large momentum transfer, 
<^ —t <^ q^. In this case, the second term of the integrand in (|8Ql) is negligible for almost 
all the range of < 1/A, and moreover, the range of z' integration is effectively limited by 
z' < (— because the modified Bessel function Kiy{\f-tz') falls off exponentially for 
z' > {-ty^/'^. In the region, < < (— t) the wavefunction of hadron shows the power 
law behavior, $(2') oc {zA)^ . Hence, (^^^^.(a/^/A) is approximately, 

gt{^t/A) ^(^] ~gt. (94) 



where g^^ is independent of t, and is of order unity. Because the factor (A/a/^)"'^" in ( 194|) is 
a rapidly changing function of iv for A^ ^ — t, this factor has an impact on the saddle point 
value V* of v integration in fl83|) . The saddle point now depends on t, as in 

(g/A, X, -t » A^) = v^^^^^f^, (95) 

In ' — 



and the saddle point approximation of the DDVCS amplitudes at leading order is given by 



r2(j^./2) 1 27rv^ y V4Vax/ V^y \^ 

Ignoring factors of order unity and the Gaussian measure of the saddle point approximation 1^ 



1 A'^"* f A\^"' ( A \ 



Im/,(x,77 = 0,t,g^)~ (-^J (^-J (^-^ ) , (97) 



,yx;j (x) \7=t) (98) 

It is customary in perturbative QCD to describe the momentum transfer t dependence of 
GPD in terms of a "form factor" F defined by 

E{x, rj, t, q^) = H{x, = 0, t = 0, g2)F(x, rj, t, q^). (99) 

Since the GPDs become PDFs H[x,q'^) for rj = and t = 0, the form factor F{x,ri,t,q^) 
should be 1 in that limit. This form factor takes account of finite size of hadrons, and is non- 
perturbative in nature. There is no way calculating the form factor in perturbative gauge 



^^The effective exponents introduced in (|89|) are still given by 7cff. = 7i/» and Aoff. = jv* — 1; the saddle 
point value of v* is now given by (^5]) . Thus, the discussion in pl^B] holds true also for this case without 
modification. 
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(a) (b) 

Figure 5: (color online) The panel (a) shows the form factor F fllOOl) as a function of — t/A^ 
in low momentum transfer {—t < A^) region. We used the scalar target hadron with A = 5, 
n = 1 explained in section |5J The three curves in the panel (a) correspond to finite x, iv* = 4, 
smaller x, iv* = 2, and small x limit, iu* = from above (yellow/light) to below (blue/dark). 
The panel (b) is a schematic picture of form factor for large momentum transfer. Dashed 
curve corresponds to a smaller x than the one for the solid curve, and short dashed curve 
to an even smaller x. The form factor eventually shows a power-law behavior for sufficiently 
large momentum transfer; this power behavior is reached for smaller momentum transfer for 
smaller x (see discussion around fllUip ). as indicated in the figure. 



theories^ and experimentally measurable form factors (such as the electromagnetic form 
factor) are sometimes used for theoretical modeling of GPD (though it should depend on 
(x, 7], q^)) and for fitting of DVCS experimental data. With the holographic set up, however, 
it can be calculated from first principle. 



Just one common form factor F is necessary for all the structure functionqffj ^i,...,5, 
because the expressions ( 188} [98]) are both not much different for Jq and Ji , at least for small 
X we have assumed so far. Taking ratio of flHTt |88| |98|) to /j(x, r/ = 0, t = 0, g^), we find that 

(100) 

The form factor largely shows power-law dependence on (— t) at large momentum transfer. 



■^"'in the approach of [57l EHl [59l [60l |6T] , power-law scaling in energy can be concluded for large (— i) region, 
but still non-perturbative wavefunctions of partons are required for quantitative results. 

^^Here, we talk of form factors describing t dependence in the structure functions, rather than the t 
dependence of GPD. Characterization of GPD in strong coupling regime is given in section where we will 
see that the form factor of structure Junctions (jlOOp can also be taken as that of GPD in the sense of (|99)) . 
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and for sufficiently large momentum transfer, 



In ( -1= W ^ In f-^ ) , (101) 



the exponent of is approximately given by —(2 A — 2 + jo)- 

The Regge behavior in the Virasoro-Schapiro amplitude (sio)"'*^°^^ does not lead to ex- 
ponential dependence on the momentum transfer t of DDVCS processes, because this fac- 
tor works as an exponential cut-off (s)^'^^/^^ for contributions from Pomeron with large 
Kaluza-Klein momentum in the direction of holographic radius. Pomeron with small u still 
contributes, but the Pomeron wavefunction '^'fj {t, z) ~ Ki^{y/^z) cuts off the IR z^l/ 
region of AdS^ for spacelike momentum transfer [3]. This exponential cut-off and the power- 
law wavefunction (jH]) of the target hadron combined results in the power-law "^^^^^ 
dependence of the form factor, as we have already seen in (IMl) ^ 

Note that this power-law dependence of the form factor is a generic consequence of asymp- 
totically conformal theories, and is independent of details of IR geometry in holographic 
models. The power is determined by the conformal dimension A of the bulk field the hadron 
belongs to. The power-law behavior in the UV conformal holographic models is due to the 
cost of squeezing the entire hadron into a size of 1/a/— t. This power-law behavior in holo- 
graphic picture, however, should not be taken as an explanation for the power-law behavior 
observed in the real world hadron-hadron elastic scattering data at large momentum trans- 
fer ^ — t, because the holographic models assuming large 't Hooft coupling even at high 
energy is not truly dual to the real world QCD. The power-law behavior expected from the 
naive power counting based on the number of valence partons [571 EH EH EOl |61] should be 
close to the truth of the power law in the <^ —t region. Despite this difference, it is an 
encouraging fact that the form factor F{x,t,q'^) of holographic models generically matches 
on to something semi-realistic (power-law behavior) at large momentum transfer A^ <^ —t, 
rather than to something totally different. There is a long tradition of using a GPD form 
factor with a power-law behavior in t at the A^ ^ —t region (e.g., [M| W7\ EB]), and it might 
be possible to provide a better theoretical foundation to such form factors. See also section 

El 

To learn more about the structure of the non-perturbative form factor ( llOOp in the strong 
coupling regime, let us examine the kinematical variable dependence of the form factor fllOOp 



■^^This mechanism is similar to the way the power-law dependence of DIS cross section was obtained for 
moderate x [2], although the integration along the holographic radius was limited by the non-normalizable 
wavefunctions of virtual photons P^t^^. (z) ([751171)) in the DIS case, not by the Pomeron wavefunction {t, z) 
as in this case. Gravitational form factor — the j = 2 moment of non-skewed GPD — also shows power-law 
dependence on the momentum transfer [H [62l [63l [64l [65] . 
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more carefully. First, the x dependence and (— t) dependence are not completely factorized. 
Although the form factor gradually approaches to power-law in with a common power 
— (2 A — 2 + jo); it does so under x-dependent condition f llOll) . Even within the region fllOip . 
the coefficient of the power-law (A/-\/^)^^~^+-'° also depends on x and g^; we find that an 
effective energy scale Aeg in F ~ [A^g/(— t)]^"^"*"-"'/^ is given by 

v^[ln(g/A)]^ 

A^jj = A X e(2A-2+Jo)lii[(9/A)/v^a;] (102) 

and in particular, A^jj decreases for small x (large W"^). Note also that this form factor and 
A^g properly take account of dependence, not just (x, t) dependence at a fixed value of q^. 

In the extremely small x limit (for a given value of g^), however, Aeff, does not become 
arbitrarily small, but approaches a finite value, A^. In this limit, the form factor greatly 
simplifies, to become 

F(x,r/ = 0,t,g2) ~ (^-^j (A^ « -t « g^), (103) 

from which x and g^ dependence has disappeared. This is when iv* = 0. 

4.1.3 Regge theory revisited 

Regge theory and the dual resonance model had a close relation in description of hadron 
scattering in late 1960's-early 1970's. Certainly the dual amplitude on 3+1 dimensions were 
not able to explain power-law behavior in the fixed angle high-energy scattering on one hand, 
and the theoretical consistency of string theory hinted a spacetime dimension higher than 
3+1 on the other. But, string theory on a warped spacetime has resurrected as a theoretical 
framework of high energy scattering of hadrons in strongly coupled gauge theories [H E], [3] . 
When the scattering amplitudes of holographic QCD are seen as an amplitude of hadrons in 
four dimensions, some part of the Regge behavior of the scattering amplitude of strings in ten 
dimensions still seem to remain. The Pomeron trajectories predicted by holographic QCD 
are not linear, on the other hand, and the complex j-plane description of hadron scattering 
amplitude of holographic QCD may be a little different from what one naively imagines from 
traditional Regge theory with some linear trajectories. The question is, then, how they are 
different. 

Preceding literatures [21 [HI [12] have already made efforts in providing complex j-plane 
description of the scattering amplitudes of holographic QCD. Our discussion in sections 
l4.1.1f[^?L2] already shows the importance of the saddle point value j = j^* of the scattering 
amplitude in j-plane in extracting kinematical variable dependence of the scattering ampli- 
tudes. It is thus worthwhile to take a moment in this section 14.1.31 and elaborate more on 
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j-plane description of the hadron scattering amplitudes in holographic QCD. It is known that 
i) holographic QCD gives rise to Kaluza-Klein towers of Pomeron trajectories, and that ii) 
they are not linear. From these properties, we will see in the following that scattering am- 
plitudes are described better by treating those trajectories individually for some kinematical 
range, while they are better described by treating a Kaluza-Klein tower as a whole for some 
other kinematical regions. Such transitions are triggered by the location of the saddle point 
in the j'-plane relatively to other singularities of the amplitude. 

Let us begin with studying behavior of the DDVCS amplitude, not just for physical region 
of {s,t) plane, but for region including real positive t. The integrand of fl7^ can be regarded 
as a holomorphic function of u and t. The denominator of the second term, /^^/(a/^/A), 
becomes zero if 



t = K\3^u,nf 



n 



1,2, 



,GN, 



(104) 



and these zeros in the denominator can be regarded as poles in the i^-plane for a given 
kinematical variable t. As we analytically continuate the integrand from real negative t -C 
— A^ in the physical kinematical region to real positive t^A^ through the upper-half complex 
plane of t, the poles in the z/-plane move and some of them show up in the region with real 
positive iv; see f ll04p . We introduce a notation tc,n = (Ajo,n)^; if ^c,m < t < tc,m+i, there are 
m poles in the lower-half z/-plane. The positions of such poles satisfying fll04p are denoted 
by i^n(t)- Thus, for A^ <^ t, the integration contour in the i^-plane needs to be deformed as 
in Figure E] (a), which can then be rearranged as in Figure 
Contributions from the poles u = z/„(t) can be written aH^'' 




E 

n=l 



where 



and 



1 + Q~^ij:.„(t) 

sin(vr>„(t)) 



r2(>„w/2) 



/3n(«p,n(t),t; 



( 



Jo + 



ap,n(i), 



= Vn(t) 



An 



77*7*Pn(cip,n(i)) lhhVnip'V,n 



(105) 



(106) 



(107) 



Although we obtained (|105p by analytical continuation of the amplitude from i < into i > through 
the upper half plane, there is another derivation. In section 13. 1[ we derived the Pomeron kernel by inserting 
the complete system ([58|) into (|54| . but we implicitly assumed that t < 0. For ^ t, however, ([58|) is not a 
complete system, as the "Schrodinger equation" — Aj(t)4'^^' = [E^ / E?)'^^^^ may have discrete spectrum as 
well [3]. Indeed, the wavefunctions (|117|) become the discrete spectrum of Aj(t) by replacing ivj with ivn{t) 
and rUj^n with a/Z, respectively. The pole contributions (|105p come from the discrete spectrum part of the 
inserted complete system. 
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dzy^Phhiz){RAye 
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(108) 

. (109) 



Under the conditions that 



j = ju (equiv. z/ = Uj) and v = z/„(t) (equiv. Vt/A = ji^^n), (HO) 

(from which j = ap^„(t) follows), u = i>j = Vn{t)i and a j dependent mass parameter 
rrijn = rrin^ is defined by 



(111) 



The factor 77*7*pn(j) can also be defined similarly to (11081 1109p by replacing Phh{z) by 
P^l?^*{z) {i = 0,1) in ( 1731 [71]). Because the wavefunction of photon depends on g^, it has 
dependence on g^. Explicitly, one can find 77*7*pn(ttp,n(^)) ~ (A/g)'^("''"(*))+^"'P'"(*\ and 

1 + e-'rji^„{t) 



n=l 



1 ~ /^\ TCop.nC*)) /-|\ap,n(i) 



sin(vrj^„(t); 

where ^„(t) = /3„(ap,„(t), t; g2)(g/A)7(«p,n)+2ap.„ jg independent of q 



(112) 
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The expression f llOSp is precisely in the form assumed in traditional Regge theory: 



A 



Regge 



j=a(i) 27rz sinvrj J - a{t) \so 



1 _|_ p-iTra{t) / 



Traditional theory was only able to assume a linear form 

a{t) = ap^o + ctp^ 



:ii3) 



(114) 



for simplicity or from fit to experimental data, but holographic QCD predicts a trajectory 
(11061) : it is approximately linear for <^ t (but not too large —t), but it is not unless 

<^ t. The residues /3n(ap,n(i), ^; Q'^) satisfy factorization condition; they factorize into 
'j^*^*Pn{(y¥,n(t)) and 7/i/iPn(«p,n(^)) holomorphically in j = an(t). This factorization is neces- 
sary for unitarity of hadronic scattering processes in 3+1 dimensions, which is not guaranteed 
a priori in a theory that is not based on a local field theory on 3+1 dimensions [2S1 EHl EOlH 
Scattering amplitudes of dual resonance model and superstring theory in 10-dimensions have 
the factorization property [TTj, and the factorization (I107P in 4-dimensions is remnant of the 
factorization in string theory on higher dimensions. 

We have so far dealt with the (analytic continuation of the) scattering amplitude /j, 
but closer connection to the classical Regge theory can be established by clarifying physical 
meaning of ap,„,(t) in (I106P and of the factorized residues /3„ oc 'jhhVnij = o:v,n(t)) (11071 
11091) . The Pomeron-hadron-hadron coupling '~fhhVn{t) (I109p looks like an overlap integration 
of three wavef unctions; two of them Phh{z) = c^Y^{z)Y are those of the target hadron, and 
the last one satisfies 



A,(t) 



t/^ + 4 



0. 



(115) 



For a special case as in (I109p . when z/, t and j are related by z/ = Vj and t = rn^^, this means 



A.Kn)--0'-2) 



a' 



[e^^^-')Ji.,(m,-„^)] =0, 



(116) 



which is the equation of motion for spin j string stat^ on the graviton trajectory in AdS^, 
when j G 2N; this momentum in the 3+1- dimensional Minkowski spacetime is on the mass 
shell t = m?- in this wavefunction, when the IR boundary condition is satisfied. Therefore, 



Factorization predicts relations among differential cross sections: dcrc\{A + B)/d(Jc\{A + C) = d(Jci{D + 
B)/d(Jd{D + C) [30]. 

The mass of the string will differ from one in the flat space. But we have shown in section [5TT] that the 
mass shift 5j can be consistently neglected within our approximation. 
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j = «p,n(^) describes the spin-mass relation of 4D hadrons corresponding to the n-th Kaluza- 
Klein mode in the spin j string state in the graviton trajectory; j = ap ,^=i(t) is the leading 
trajectory, and j = a^^nit) with 2 < n < m become daughter trajectories [3]. Normalizable 
wavefunctions of all the 4D hadrons in the trajectory j = (11 101) are reproduced from 

a single wavefunction 



defined for j G C; the wavefunction for j G 2N is simply the special case of the one above. 
One might refer to (I117P as a wavefunction of the n-th Pomeron trajectory. The Pomeron- 
hadron-hadron coupling '~fhhPn{j) is given by an overlap integral of two $'s along with this 
wavefunction ijjn\z) of the n-th trajectory. That is, 



When t < on the other hand, all the poles in the lower half z/-plane have moved back 
to the upper half plane, and the scattering amplitude li is given by a continuous integration 
of u real axis. The u integration in (155]) can be converted into j integratioE0 through a 
change of variables j = j^, but there is no pole in this j integration, unlike in the traditional 
form of Regge theory amplitude (11131) : the DDVCS amplitude (175 | 15 ^ 1571) that we studied in 
sections l4.1.1f[^?L2]fo r t < do not seem to be based on the Pomeron pole exchange idea, 
at least apparently]^ If so, while the coefficient of the (iy2y=a„(t) f^c^r for tc,n < t, 1hhF„{j), 
can be characterized as hadron-hadron- [Pomeron j = an{t)] three point coupling (IllSp . how 
should we characterize the coefficient of the (VT^)-''^* factor for t < tc,i, gi^*{\^^/ A), in the 
absence of Pomeron exchange picture? 



Field theory dual language is useful in characterizing the factors Cj* , A^h and g^^,{y/^/ K) 
in dTHHHO])- Let us define 



for a parameter e that has a dimension of [length]. Then one can see that the two factors 
[CW(j, q)]i/e and [Ahh{i, W"^, t)]i/€ correspond to OPE coefficient of a "twist-2" spin j operator 
renormalized at fi = 1/e, and matrix element of the spin j operator renormalized at 

^°The integration contour becomes the one in Figure [7] (a). 

It should also be kept in mind that the sublcading contribution comes from the exchange of the second 
Pomeron trajectory (the n = 2 term in (|105p for <c.2 < t, which is power suppressed, x (A^/M^^)"!^*'""^'*^ 
Corrections to (|87l [96]) on the other hand, are suppressed only by ^/X / ln{W'^ / A^) . 




.2 1 1/2 



(117) 




(118) 
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[72| [73| El] . To be more precise, we can define 



[^hh¥* {j,t)] 



t/A) X (eA) 



2A 



(6A) 



7(j) 



t/A) 



t/A) 



Iiu,iV-tz) 



(120) 



by pulling out (VT^)-'' ~ {'^q-vY from [^/^^(j, W^^, t)]i/^. It is easy to see that the short distance 
scale R drops out from (11201) . and [r/j/ip*(j, t)]^i^ can be expressed only in terms of low-energy 
data and renormalization scale = 1/e. Now one can see that 

[C«(j, g)]v. [A.(J, W\ t)] ~ (^ (^(^ ■ ^)'' t^/^'^^-O-' 



je2N 



(eg)^O) (g^ 



2^,. [^M.-'-^M.] (MP2)|[r'^^-'^^]i/JMpi))- 

(121) 



Therefore, \^hh¥*U^'t)]i/e regarded as the coefficient of [p^^ ■ ■ -p^^] of the matrix element 
of the target hadron h with insertion of a twist-2 spin j operator \T^^"'^^]ii^ (that is, spin j 
form factor or reduced matrix element). Thhf*{j^t) is now defined for arbitrary j G C in the 
expression above, not just for j G 2N. 

Using Kneser-Sommerfeld expansion of modified Bessel functions]^ one can see for arbi- 
trary complex j and t that 



[rfthp*(j,t)]i/, = Y 



IhhWnij) [^) [^) 



7(i) 



t-m2„ AJ-2 

n=l J'" 



1/2 



(123) 



The spin j form factor Thh¥*{j,t) to be used in describing the DDVCS amplitude for t < tc,i 
and the hadron-hadron-[Pomeron j = an{t)] three point couplings ■yhh¥n{j) are related as 
above. Such a relation between a form factor F{t) and three point couplings of Kaluza-Klein 
hadrons gthn, 

oo ^ 
^it) = y , jghhnFn, 



n=l " 



For < w < < 1 and arbitrary complex numbers /i and r, 



TT Jf_,{Tw) 

4 J^(t) 



1 



)l=l 



(j/i,ri)^ 



(124) 



(122) 



This is regarded as the Green function on AdS^ with an infrared cut-off, seen as a propagation of a 5D field 
(left hand side), or of a Kaluza-Klein tower of 4D fields (right hand side). 
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has been known for conserved currents [HOE]. Here, f ll23p establishes such a relation 
simultaneously for arbitrary j G C, not just for a given spin, and that makes it possible to 
"sum up" Pomeron exchange amplitudes labeled by the Kaluza-Klein excitation level n. 

* * * 

Now we know that the amplitude is given (approximately) by Kaluza-Klein tower of 
Pomeron pole exchangj^ for large positive t <^ t), while the entire amplitude as a whole 
is approximated by the saddle point method for large negative t {t < tc,i)- Since tc,i = (Ajo,i)^ 
in the hard wall model is positive, the latter should be applied for all the physical region 
t < 0. Such subtle things, however, depend on details of infrared geometry of holographic 
models, and cannot be regarded as a robust prediction of holographic QCD. 

Imagine, for example, a holographic model corresponding to an asymptotic free running 
coupling constant (e.g. [75]). With the AdS^ curvature changing over the holographic radius, 
the entire spectrum of the Schrodinger equation — Aj(t)\E'£^(t, z) = {E^/ R^)'^^J [t, z) becomes 
discrete [3]o It is now more convenient to write the amplitude in complex j-plane, rather 
than in z/-plane. The branch cut at j < jo in the hard wall model is replaced by a densely 
packed poles along the real j axis, and there is no special value like tc,i for some critical change 
in the spectrum. The spectrum may change for different values of t, like in Figure [7] (c) and 
[7| (d, e), but they are not different qualitatively. Thus, there is nothing wrong a priori in 
choosing the integration contour as in (d) or (e), even for a given spectrum. 

There exists a convenient choice of the contour in the j-plane, however. That is to let the 
contour to pass the saddle point j = j^* , and treat all the poles to the right of the saddle 
point (i.e., j^* < Re j) as isolated discrete contributions. All the rest are treated as if they 
formed a continuous spectrum as in the hard wall model. As we will see in the following, 
whether there are some poles remaining to the right of the saddle point (e.g.. Figure [7] (b, d, 
e)) or not still works as a criterion for various physical transitions, independently of detailed 
difference in various holographic models. 

As we are talking about closed string amplitude, contributions from a single trajectory cannot be regarded 
as their ^-channel exchange; they can be regarded as a sum of t-channel and u-channel exchange, however. 
See the review material in section [3. II 

Although we talk of "asymptotic free" running, it is safer to consider gravity dual models corresponding 
to gauge theories with asymptotic free running only up to some energy scale, above which they become 
strongly coupled conformal theories (e.g. [7S|). In this case, although there are many bound states in the 
spectrum of the Schrodinger equation, the continuous spectrum still remains. Thus, the densely packed poles 
are still followed by a branch cut along the real j axis in the small Rej region. As long as one pays attention 
to a certain region in the j-plane, rather than to formal difference in the literally large negative Rej limit, the 
such models pass for "asymptotic free" gravity dual models. Whenever we refer to asymptotic free gravity 
dual models in this article, this must be understood. 
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(a) 


(b) 






— 










• • — • — o • 



(c) (d) (e) 



Figure 7: Singularities and integration contours in complex j-plane. Hard-wall model is 
assumed for (a) t <^ tc,i and (b) tc,i ^ while a holographic model for asymptotic free 
running coupling is assumed (c) with a smaller t and (d, e) with a larger t. The panels (d) 
and (e) differ only in the choice of contour. Black dots are poles, wiggling lines in (a, b) are 
branch cuts, and open circles in (a-e) denote saddle points of the amplitude on the complex 
J -plane. 



38 



The saddle point in the j-plane is understood as a consequence of two competing effects 
in the DDVCS. One is the (H^^/A^)-^ ^ {q^/xK^Y factor, which is large for large real j. 
The coupling of spin j string state with a virtual photon, on the other hand, behaves as 
{A/qy^^"'^^\ where 7(j') is the anomalous dimension of the operator corresponding to the 
string state. This second factor is small for large real j. These two factors combined, ( 19T|) . 
forms a saddle point j* in the j'-plane for W ^ q ^ A (meaning x ^ 1); it is quite generic 
for holographic models with <^ a' that the large j behavioi@ of the anomalous dimension 
is 7(j) ~ [-\/Aj]^/^. (We neglect here correction from the factor of power of y/t/A.) Thus, 
in the right-hand side of the saddle point j* in the j-plane, j* < Rej, the larger the real 
part of a pole, the larger the contribution of the pole is. Writing down the amplitude as a 
sum of these poles, starting from the one with the largest real part of Re j to the ones with 
smaller Re j, we obtain a finite term sum that gives a good approximation to the amplitude. 
Individual contributions from the poles in the left hand side of the saddle point j*, on the 
other hand, becomes larger and larger as the real part of j of the poles become smaller. Thus, 
their sum does not make sensf^; all contribution in the left hand side should be treated as 
an integration on a contour as in Figure [7] (d). All of these contributions combined can be 
evaluated by the saddle point method, as we have presented by using the hard wall model, 
and their contribution is 



This contribution is even smaller than the one from a pole whose real part is even slightly 
larger than j*. This is why it is convenient to take the contour in j-plane so that it passes 
the saddle point as shown in Figure [7] (d). 

Then we can understand that there is a transition depending on whether the saddle point 
value j*{x, q, t) has a larger real part than that of the leading singularity (one in the j-plane 
with the largest real part), j = api(t). It is easy to see this in Aefr. (which we have already 
discussed in section 13.1.11) . 



*^We are not talking about j as large as \/X, however. 

In "asymptotic free" gravity dual models in the sense of footnote the series of poles stop at a certain 
small value of Re j, and a branch cut starts toward large negative Re j. Thus, the contour can be chosen so 
that the amplitude is given by a sum of all the individual pole contributions and an integral around the branch 
cut. The sum is thus, formally, well-defined in this sense. It will be obvious, because of the discussion in 
the main text, however, that the resulting many-term summation contains large cancellation between the cut 
contribution and the pole contributions, and is not as practically useful an expression as what we described 
in the main text. 




(125) 




q < qc{x,t), 
q > qcix,t), 



(126) 
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a2 q3.{x,t) 




(b) 



Figure 8: Transitions (crossover, to be more precise) in (a) Acff, as a function of for fixed 
{x,t), and (b) in [d[lm Ii\'^/dt\ as a function of t for fixed 



where the transition is induced (assuming that 7(j) is a decreasing function of j along the 
real axis) at 



dj 



dj 



ln(l/v^x) 
ln(ge/A) • ^'''^ 



j=ar.i{t) J j=j*{q=qc) 

Schematically, it behaves as in Figure |8] (a). In the following, we will refer to the two phasej^l 
as 

• Low-KK Spin-sum phase or Leading pole phase (leading singularity phase), where 
the leading singularity has a larger real part than the saddle point, and 

• KK-sum Spin-sum phase or Saddle point phase, where the saddle point value is 
larger than the real part of the leading singularity in the j-plane. 

It is intuitively obvious in gravity dual descriptions that the saddle point phase is realized 
(and such observables as 7ofr. and Aofr. are controlled by the saddle point value j*) for larger 
g^, not in the other way around; the holographic wavefunction of photon with larger virtuality 
makes the photon- [spin- j string] couplings weaker, and the couplings for larger spin states are 
affected more severely than those for smaller spin states, because of the stronger power-law 
behavior of the holographic wavefunctions of the higher spin states. See (179|) and definition 
of 77*7*pn(j)- The factor (A/g)^-'^'^^-^^ comes from there, and consequently, the saddle point 
value j* shifts to the right in the j'-plane to enter into the saddle point phase. See Figure O 



'^''The KK-sum Spin-sum phase (saddle point phase) is divided into two phases, when the real part of the 
scattering amplitude is studied. One of them is still called KK-sum Spin-sum phase (saddle point phase), 
and the other as KK-sum low-spin phase (spin-2 phase). See [11] and section l473l of this article. 
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ln(_q/A) 



Figure 9: Phase diagram (at a given value of t) in the strong couphng regime. I: leading 
singularity phase, II: saddle point phase with 7eff, < and III: saddle point phase (spin-2 
phase) with < 'jes.- A region where higher genus amplitudes are just as important as the 
sphere amplitude, phase IV, is characterized by the condition |;y'^p'i™'^(6 ~ 1/^)1 ~ C'(l), just 
as in [11] . The leading singularity phase I is absent in the forward scattering t = in the hard 
wall model, as in [11]. This phase diagram in the strong coupling regime looks quite similar 
to that of QCD. The boundary lines between the phase I and II, and II and III, however, are 
nearly vertical in the strong coupling regime, rather than being nearly horizontal as in the 
weak coupling regime. 

In fact, the transition between the two phases is not a singular phase transition but a 
crossover, unless we literally take a small x limit. As long as \n{l/x)/\/X remains finite, 
the saddle point approximation is never exact, and the notion of the saddle point itself 
should be accompanied by a width proportional to [ln(l/x)/V^]~^/^. Although we defined 
the two phases by simply comparing the real part of the leading singularity j = ap,„=i(t) 
and the saddle point j*, we should also keep in mind that the n = 2 and higher Kaluza- 
Klein contributions give rise to finite corrections (for finite \n(l/x)/\/X) to the leading pole 
contribution in the leading pole phase; the corrections (and cancellation) may be sizable for 
negative t and finite \n{l/x)/\/X. For those reasons, the transition between the two phases 
can be a singular genuine phase transition only in the ln(l/x)/A/A — ?■ oo limit. 

Phases are determined by kinematical variables {x, t, g^), and kinematical variable depen- 
dence of all the observables, not just Acff., will be different for different phases. Real part 
to imaginary part ratio of the scattering amplitude, which we study in section 14. 3^ is an 
example. Note also that the saddle point value j* depends on momentum transfer t, not just 
on (x,g^), and the leading Pomeron pole also depends on t. Thus, the crossover should also 
be induced by the kinematical variable t (at least somewhere in the t-plane). In section 14.1. 4[ 
we study another observable, t-slope parameter B{x,ri = 0,t,q^), and discuss its crossover 
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behavior (like in Figure |H] (b)). 



4.1.4 Slope parameter of the forward peak 

We have so far focused on various features of the DDVCS amphtude that are robust and 
do not depend on detailed difference of holographic models. Individual holographic models, 
however, have full control over non-perturbative aspects of hadron physics, and more (possibly 
model dependent) information can be extracted. The expression ( 187|) . for example, tells us 
how to calculate full t-dependence of the (imaginary) part of the DDVCS amplitude, not just 
for t = and the <^ —t asymptotic region. 

The slope parameter of forward peak (t-slope parameter) in the elastic scattering of two 
hadrons, 



d , \daci 



s,t) 



:i28) 



dt 

is an observable that characterizes the transition from t ^ region to < -t region, and 
has been measured for p+p and p+p scattering at various energy scales (e.g., [30]). The slope 
parameter has also been measured in HERA experiment for 7* +p — 7+p scattering (DVCS) 
[76| U7\ [78] . In this section I4.1.4[ we will first use the hard wall model to derive an explicit 
prediction of the slope parameter in the photon-hadron scattering at vanishing skewedness 
1] = (DDVCS process) for simplicity. Later on, we will discuss how much predictions on the 
slope parameter could be different for different holographic models, and discuss the possible 
crossover to be seen in the slope parameter (for some holographic models). 
We define the slope parameter in non-skewed DDVCS b} 
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Bi{x,7] = 0,t,q^) = 2 ^^Inlm Ii{x, 7] = 0,t,q^). (129) 

Let us first work out the prediction of the slope parameter in the saddle point phase, and 
then present the result in the leading pole phase later. 

In the saddle point phase (KK-sum spin-sum phase), the factor (l/xY* x (A/q)'^^^*^ in 
(1871) does not contribute to the slope parameter, because the saddle point value j* = jy* does 
not depend on t for < (— t) < A^. The slope parameter comes entirely from Pomeron- 
hadron form factor \^hh¥*U* ^ fl'iv* ("\/~^/^)) which is regarded as a "spin j = j* form 
factor" . The Pomeron-hadron-hadron coupling behaving like a form factor is similar to the 



We could define the slope parameter B by the absolute value not by the imaginary part Im/i as in 
()129p . This choice, however, makes no difference as long as {—t) < A^, and 7eff. < (the condition for the 
real part to be outside spin-2 phase); see section 4.3.1. This is because in this region the real-to- imaginary 
ratio of li is not dependent on t, p63p . It is very likely that the most of the kinematical reach of DVCS 
measurement in HERA [76l [77l [78| is also in this region. 
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Figure 10: (color online) Slope B of the forward peak in DDVCS. The dimensionless value 
5 X is shown as a function of iu* f l86p : from top to bottom, blue (solid line) curve is 
for v^^/A ~ 0.01-0.1, red (long dashed) one for v^^/A = 1., yellow (dashed) one for 
\/^/A = 3., and green (short dashed) one for \/^/A = 6. We used the wavefunction of the 
target hadron ([8]) for the first excited (n = 1) mode with conformal dimension A = 5, and 
the 't Hooft coupling was set to \/A = 10 in this calculation, t-dependence is very weak for 
< -t < A2. 



idea advocated in [791 EQ], but the form factor turns out not to be precisely the same as the 
electromagnetic (spin 1) or gravitational (spin 2) one. There exists a notion of form factor 
^hhF*{j, t) that is holomorphic in spin j, and the one with the saddle point value of j = j* is 
relevant for the DDVCS amplitude; the saddle point value j* = jy* is determined by (x, t, q^) 
as we have already seen in sections I4.1.1H¥X^ 

It is straightforward to calculate the slope parameter in the saddle point phase, 

d 



B{x,r] = 0,t,q' 



2^ In [ThhF^ 



(130) 



by using the explicit expression of the form factocj (11201) . The result is shown in Figure [TOl 
The larger the spin j* (and hence iu*), the smaller the slope. At t = 0, the slope parameter 
is now the same as the "charge radius square" of the hadron under a "spin-j* probe". 

When the slope parameter is seen as a function of and W"^, it changes only through 
the change in the saddle point value j* = jy*. Because 



>0, 



<0, 



(131) 



91n(gVA2) ln(gVA2) ' d\n{W^ / h?) ln(l^VA2) 

the slope parameter decreases for larger \a{q^ / h?), and increases for larger \n{W'^ / h?) in this 
saddle point phase prediction. It must be more sensitive to \n{q^ / h?) than to \\i{W'^ / K^) for 



A grossly incomplete list of literatures on fixed spin form factor in holographic methods will include 



43 



small X [ln(l/x) ^ ln(g^/A^)], because of the denominators in f ll3ip . That is, we find the 
following qualitative prediction of the saddle point phase@ 



dB 



<0, 



OB 



> 0, 



dB 



d\n{W/A) 



dB 



d\n{q/A) 



(132) 



d\n{q/A) ' d\n{W/A) 

There are three remarks here, before we move on to discuss predictions on the slope 
parameter from the leading pole phase. The first remark is on the size of the slope parameter 
of the forward peak. It is of order in this prediction of the saddle point phase, and is 
not tied to the slope of Pomeron trajectories; the asymptotic slope of the trajectories is of 
order (l/A^) x (l/VX) (see f lTSD l 

Secondly, let us discuss the momentum transfer dependence of the slope parameter within 
the saddle point phase. All the physical kinematical region t < is in the saddle point phase 
in the hard wall model, and the phase appears at least for sufficiently negative t {t < tc,i) 
for any other holographic models. Because of the power-law behavior of the form factor in 
UV conformal theories, the slope parameter behaves as oc l/(— t) for large (— t). The slope 
parameter does not diverge toward t — 0~ (if all the range of t < is in the saddle point 
phase), however, and it approaches a finite plateau value of order 1/A^ instead. Whether 
there is a range of t (in small \t\) where the amplitude shows exponential fall-off in t is about 
the stability of the plateau value of B for a finite range of t, and about the plateau range and 
plateau value. This is a quantitative question whose answer depends on details of holographic 
models, and we do not discuss more about this question in this article. 

Finally, it is important to note that the argument above on the W"^ and dependence of 
the slope parameter B fll32p can be applied to other observables of the DDVCS amplitudes, 
as long as they depend on and W"^ only through the saddle point value j* . Their ln(g^/A^) 
and ln(iy^/A^) dependence are always opposite, and the ln(g^/A^) dependence must be 
stronger at sufficiently small x, because the saddle point j* has three properties: 



>o, 



dr 



<0, 



df 



d\n{WlA) 



dr 



91n(g/A) 



(133) 



ain(g/A) ' d\n{W/A) 

This argument can be applied, for example, to Aefr. and 7eff.. 

Let us now move on to the leading pole phase. Although all the physical kinematics 
t < is in the saddle point phase in the hard wall model, there may also be a range of the 
leading pole phase within the physical region in some holographic models, such as asymptotic 
conformal models with negative tc,i and asymptotic free models, as we have seen in section 
14.1.31 In such holographic models, the crossover between the two phases may be observed in 



^°In this article, we assume the generalized Bjorken regime, ([51 exponentially small x (|53p . and large 
't Hooft coupling, with an extra constraint j* < 0(1) (i.e., iv* < \^/^). 



44 



the slope parameter within the physical kinematic range t < 0. The slope parameter in the 
leading pole phase comes mainly from the factor 

A2 



77*7*Pi (i) 



i=op,i(t) 



1/ 



j=ap,i{t) 



(134) 

in the small x and large region; the t-dependence in the Pomeron coupling 77*7*Pi(ap,i(t)) 
is now ignored. Thus, we have 



B{x,r] = 0,t,q^ 



' dt 



ln(l/x) + ^ 
d] 



ln(A/g) 



i=«p,i(t) 



For sufficiently small x, this is of the order of 



B 



1 



(135) 



(136) 



■ 

This is larger than the slope parameter in the saddle point phase, B ~ (^(l/A^), in the small 
X regime fl53|) we have been studying in this article. Thus, the crossover between the two 
different phases are induced schematically as in Figure [H] (b), where the large slope fall-off 
of the scattering amplitude in larger (more positive) t breaks into smaller slope behavior in 
smaller (more negative) t, when [c^7/c^j] |j=Qp comes close to [d-^ /di]\j=j*. 

Even in asymptotic conformal theories with negative tc,i and in asymptotic free theories, 
where the leading pole phase may exist in the physical kinematical region t < 0, one always 
enters into the saddle point phase for sufficiently high g^; this behavior is understandable 
just like in the case of Aefr.. 

It is an interesting question, at least from a theoretical perspective^ and also from the 
context of fitting experimental data, whether such a crossover should be expected within the 
physical kinematical range t < 0. The answer is "no" in the hard wall model. But, the 
answer depends on infrared geometry of holographic models, and a robust conclusion cannot 
be drawn only from the experience in the hard wall model. It is interesting to know the 
answer to this question in holographic models whose infrared geometry is fully faithful to the 
equation of motions of Type IIB string theory [8] 



4.2 Impact Parameter Dependence of the Imaginary Part 

One can take a Fourier transform of Im Vi{x, rj, t, q'^)'s with 77 = 0, to obtain distributions in 
the impact parameter space. Such distributions are interesting on their own, because they 

^^In reality, unitarity limit is reached and l/iVc-suppressed contributions (that we ignored throughout in 
this article) also become important at sufficiently high energy. Such high energy region is described as phase 
IV in figure [9] 
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show transverse spacial distribution of longitudinal momenta in the target hadron ^T\. Also, 
unitarization of high energy elastic scattering amplitude and DIS cross section needs to be 
discussed in the impact parameter space, because all the partial wave amplitudes, i.e., for any 
values of impact parameter, should be unitary (e.g. |30]). Although the impact parameter 
dependent amplitude of hadron-virtual photon elastic scattering has already been formulated 
in [H [ini [m [I2l [H] in holographic calculations, we elaborate more on this in the following, 
and find a result fll51H153"]) and the phase diagram in Figure [T31 which certainly refines the 
results that are already found in the literature. 

The impact parameter dependence, and hence the momentum transfer dependence, comes 
mainly from the Pomeron kernel. The profile from the Pomeron kernel is extended over the 
size of order bA ~ [ln(l/a;)/V^]^^^ ^ 1, as we will also see explicitly later. Although the 
virtual photon wavefunction also have t = — dependence through = Q*^ + ^^/4, this 
dependence is relevant to the transverse profile only in the short distance of order 6 ~ 1/q, 
so we neglect this contribution. Thus, the impact parameter dependence purely comes from 
the Pomeron kernel. Fourier transform of the kernel is 

/C(s, b, z, z') = / — — e-^-''/C(s, t = -A\ z, z') = — t/AVo(6A)/C(s, t, z, z'). (137) 
J (27r)^ An Jo 

The confinement effects (finite A effects) are crucial for the impact parameter profile at 
long distance, bA ^ 1. The profile in hard wall model shows quite different behaviors from 
one in conformal theory (corresponding to A — 0) [12]. Because we are interested in hadrons 
in confinement, we examine the explicit form of the Pomeron kernel including the second 
term of f lTB]) . in detail. 

As we have already mentioned in section 3 and section 4.1, the integrand of the Pomeron 
kernel [7H]l F^ is holomorphic (except some singularities) not just in spin j and anomalous 
dimension {iu — j), but also in momentum transfer t, except at poles t = (Ajjj,„)^ = {rrin^y 
in fll04p . Thus, one can rewrite the kernel for t = — A^ < on the real negative axis as 



K:(s,t,z,z') = ^ [ 
2m Jc 



^^,/C(M^i^ (138) 

t'-t ^ ^ 



using the function /C(s, t', z, z') holomorphic in t' (except the poles) and an integration contour 
C which goes around the non-positive real axis counterclockwise (Figure [TT] (a)). Its impact 



^^The "imaginary" part means the imaginary part when both s and t are real. The "imaginary part" and 
"real part" of the kernel separately become holomorphic functions of complex valued t. 
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■ 



A'(i^.,2 



(a) (b) 

Figure 11: (a) contour in the t-plane before deformation, along with location of singularities 
(for some fixed v), and (b) the contour after deformation. 



parameter space representation becomej^ 



]C{s,b, z, z') = / dt 

2m Jc 



/^l/2^ 



27r 



-/C(s, t', z, z'). 



(140) 



Because of the holomorphicity in the complex t'-plane, the contour C can be deformed as 
in Figure [TT] (b), picking residues at the i^-dependent poles in the t'-plane; we can do this, 
because the integrand (for given v) vanishes exponentially at \t'\ — > oo. One can show that 



/C(s, 6, z, z) 



+00 



dv iv 



sin(7rj,) J r2(jV2) V 4 



The imaginary part [resp. real part] of IC{s,b, z, z') is obtained by taking imaginary part 
[resp. real part] of [1 + e"'^*-'"], as the Fourier transform works separately for the imaginary 
part and real part in 01371) . The n-th term corresponds to the n-th term of fll05p . 

Let us examine the profile of the imaginary parts of /o,i in the impact parameter space. 



The following relation is used: 

dA 



- 2Jo(6A) 







-TT < arfft' < TT. 



(139) 
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We will focus on a region 



6>A"\ (142) 

where the series of contributions from Kaluza-Klein tower of Pomeron trajectories fll4ip 
converges quickly; since the real part of m^f ■* is of the same order or larger than the hadronic 
scale A at least for moderate value of G M and n G N, we can approximate 



Koimi%) ^ J-^e-"^-'". (143) 
V 2mrb 

Therefore, in this region, f ll42p . the amplitudes fll4ip are dominated by the leading (lowest 
Kaluza-Klein) Pomeron trajectory, n = 1, and the contributions from the higher Kaluza- 
Klein trajectories are suppressed by of order e"^^^^'-'^/^^". 

At small X, (ISTl 1531) . the impact-parameter dependent Pomeron kernel can also be evalu- 
ated by the saddle point method. In addition to the (a's)-'" and KQ{bm^^^) factors, another 
factor Ji^{m^i^ z) may also give rise to a large dependence in u; dominant contribution to 
scattering amplitudes come only from small z ~ 1/g region in DDVCS as well as in DIS, and 
hence Jiy{rn"^\) ~ {m^^^ . The saddle point z/*(g/A,x;6) is, therefore, determined by 



zu*{q/A,x;b) = ■ (144) 

i In ( ^ 1 



Therefore, we find that 

lml,{x,v = 0,b,q^) [tT.) [l ) ^ 



Im Ii{x,r] = 0,t = 0,g2) ^^y.*(,t=o) f^y-''(t=o) 



(145) 



Note that the saddle point value iv*{b) in the numerator is determined by (I144p . while iv*{t = 
0) in the denominator by fl86|) . The x-dependent transverse profile of parton distributioij^ 
is given by Im /^(x, r] = 0,b, q^) dxd?b. 

The transverse profile above is not a simple function, especially because rrL{ depends 
on X, q^ and b through (11440 . In order to extract the qualitative feature of the profile, let us 
first consider situations where the saddle point value iv*{b) is much smaller or much larger 



The impact parameter dependent Pomeron kernel can be used also to determine the phase shift of 
the scattering amplitude in the impact parameter space, x('S, 2,2') [lOj . The relation between them is 
Ii{x,ri = 0,6,(j^) ~ (sN^) x('Si ^')U~i/g,2'~i/A- section l4?3l where we discuss the real part of the 
scattering amplitude, we will study the 5-dependent phase shift, rather than 6-dependent 
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than unity. Then the following expansion of the zeros of Bessel functions of order iu, jiy^ni can 
be exploited in order to understand how the saddle point value iz/*(g/A,x; b) is determined 
by x, and h through fll44p : 



iiu,n - jo,n + Cnif + 0{{ivf), 

55| 



(146) 
(147) 



with positive numbers c„ and dn of order unity 
when 



ln(g/A) 



ln(l/xv^)/VA 



^ 1 and 



The saddle point stays within |z/*| ^ 1, 
6A 



ln(l/xyA)/v^ 



< 1 



(148) 



and the saddle point is 



ih'*{q/A, x; h) 



ln(g/A) + cihK 



(149) 



On the other hand, when either of the left hand sides of fll48p is much larger than unity, the 
saddle point becomes |z/*| 1 and 



ih'*{q/A, X] h) 



ln(g/A) + hK 
V\ \Vxx) 



q/A_ 
\/\x 



(150) 



Now it is easy to see the transverse profile for large q^ satisfying In (;|) ^ In 
this situation corresponds to (x, q^) where 7efr. is positive, that is, PDF decreases in DGLAP 
evolution, as we saw in section 14. 1.11 



In 



Im /j(x, r] = 0,b, q^) 
Im /j(x, ?7 = 0, t = 0, q^) 



A^ exp 



&Aln(|)+l(6A)^ 



-6A 



[6A«ln(i)] 



e [In (I) < (6A)] . 

(151) 

The profile remains linear exponential in the impact parameter h for smaller 6, until it turns 
into Gaussian for larger h. 



Let us now study the transverse profile for smaller q^ satisfying In (-1) <^ In 
instead. This condition on (x, g^) corresponds to iv*{^ = 0) ^ 1, when the PDF sti 
under the DGLAP evolution, 7eff. < 0. The saddle point value iu*{q/ K,x]h) becomes large 



g/A 
\/\x 

1 increases 



55^ 



7r/2 for large n. 
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for large b which violates the second condition of ( 11481) . In this case, by using fllSOp . the 
transverse profile turns out to be Gaussian approximately: 



Im r] = 0,b, q^) 
Im Jj(x, ?7 = 0, t = 0, g2 



9/A 



for — = In 



g/A 
a/Ax 



< (6A) 



(152) 



In a range of smaller impact parameter where iv*{h) < 1, however, the expression cannot be 
simpler than 



Im Ii{x, 7] = 0,b, q^) 
Im /j(x, 7] = 0, t = 0, 



3m 



A^e- 



for (6A) < — In 
V A 



g/A 
-\/Ax 



53) 

The transverse profile in this range of impact parameter h is approximately linear exponential 



in h (c.f. |T7]), with the mass parameter rn( 



iu'{q/A,x;b)) 



changing slowly in {q/A,x,b). This 



approximate linear exponential profile fll53p smoothly turns into the Gaussian profile 01521) . 
because of the b dependence of the saddle point iiy*{q/ A,x;b). All of these resultqfj are 
summarized in Figure [13] (a). 

Regardless of whether is large (11511) or small (11521) . the imaginary part of the (sphere 
level) DDVCS amplitude shows Gaussian profile for sufficiently large impact parameter b. The 
Gaussian profile in the transverse direction has been used for phenomenological fit of GPD, 
and the holographic calculation above justifies the Gaussian ansatz at least for sufficiently 
large impact parameter. 

This agreement between the holographic calculation and the conventional phenomenolog- 
ical ansatz is not an accident. The phenomenological Gaussian ansatz is a direct consequence 
of a linear trajectory (I114p in the traditional Regge ansatz (e.g. [30] )• Although the Pomeron 
trajectories (11 101) in holographic QCD are not linear at all, they become approximately linear 
at large t and j (but not too large j [S2]), because of (I147p : 



1 



(154) 



Since the saddle point for large b is determined by the large iu behavior (and hence by the large 
t (and large j) behavior) of the Pomeron kernel, the asymptotically linear trajectory (I154p 
gives rise to the same large-6 behavior (that is, Gaussian profile) as in the phenomenological 



Note that for (6A) ^ ln{l / {^/Xx)) , neither (|15ip nor (|152|) are rehable, as we have used the form of 
jrii^) that is rehable only in ^ ^/X. For such a large impact parameter, the saddle point value of i' is not 
within this range. Similarly, for ln(g/A) :$> / (^/Xx)) , (jl5ip cannot be trusted. 
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ansatz with a linear trajectory. The width-square of the Gaussian profile should be given by 
using the asymptotic slope parameter of the Pomeron trajectory j = ap,i(t) in (1154^ : 

which is indeed the case in (11511 1152p . 

The Pomeron trajectories in holographic QCD are far from being linear, however, except 
for the large j-t asymptotic region. As a result, the transverse profile may well be different 
from the Gaussian profile at smaller b (the region where bA <^ \n{l / x) / ^/X) . Indeed, in this 
region, the results (1151^ I153p show that the profile is not Gaussian, but approximately linear 
exponential in the hard wall model. Interestingly, the mass scale of linear exponential profile 
for small in (I153p . rrin depends on kinematical variables, and is different from 

mass eigenvalues of any 4D hadron states; this linear exponential profile is not associated with 
a t-channel exchange of a single particle, but emerges after summing up all the stringy states 
with different spins in the n-th Kaluza-Klein trajectory. The mass scale 'm^n=i^ appearing 
in the small x limit (for fixed and b) is smaller than any one of the mass eigenvalues of 
the stringy (and graviton) states in the trajectory, and the range is longer than a simple 
t-channel exchange of a glueball. 

Although the expressions (I151H153"|) rely on the hard wall model, most of its qualitative 
aspects are expected be in common with other holographic models. The transition from (I152[ 
I153P to (I15ip for larger g^, for example, will be induced in various models. The saddle point 
is dragged to larger j, iu and t for a given value of b, just like we discussed in section 14.1.31 

4.3 Real Part of the Amplitudes 

The imaginary parts of the structure functions are related to the GPDs, whose Fourier 
transforms are interpreted as the distributions of partons in the transverse spacial directions. 
Although the real part of the (D)DVCS amplitude does not have such an interpretation, it 
still contributes to the DVCS cross section. 

The large impact parameter b behavior of the real part has been discussed in [12], and 
the behavior in the small b region (but not as small as bA < 1) is covered by [T7]. Overall 
normalization is found in |TT] . The following discussion provides a clear description of how 
the small b behavior of {^7\) smoothly turns into the large b behavior of [12], as well as 
careful interpretation of the physics behind this phenomenon. This subject has, in fact, quite 
a long history (e.g., [69]), but we hope that the following presentation using the saddle point 
value j* as a key concept helps clarify things a little bit. 
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4.3.1 Momentum Transfer Dependence 

It is convenient (and customary) to use a variable 

i. „ 2N _ Re /i(a;,r/ = 0,t,g2) 

Im /i(a;,?7 = 0,t,g^) 

in characterizing the real part of the amplitude, as we already know the behavior of the 
imaginary part. The value of p{t) at t = is often denoted by p in elastic scattering of 
hadronsO The real part of the amplitude is obtained by simply taking the real part of 
[1 + e"'^*-^] in the Pomeron kernel (1641) . The ratio p{t) is simply given by 

p(t, X, r/ = 0, q^) = tan (| (j*(x, r/ = 0, t, q^) - l)) , (157) 

using the saddle point value j* = iu*[x,q^,t) in the hard wall model [TT] for the entire range 
of physical momentum transfer t < 0. The j-v integration in the kernel can be evaluated by 
the saddle point method for small x (1531) . unless q^ is too large to satisfy 

iv* < 2 (equivalents j* < 2). (158) 

The expression for p{t) above remains valid in other holographic models, as long as (x, t, g^) 
is in the saddle point phase. p(t) ^ 1 is predicted in gravitational dual in general, because 
j* is closer to 2 than to 1 for A > 1 [Tn[T2]. 

In holographic models other than the hard wall model, there may be a leading pole phase 
in the kinematical variable space (x,t, g^). There, the ratio p{t) is given by the expression 
f ll57p with the saddle point value j* replaced by the leading Pomeron pole ap^i(t) [30] . 

The relation f ll57p follows immediately from derivative analyticity relation [SH], if j* — 1 
is understood as X^s., the effective exponent of s = oc 1/x of photon-hadron scattering; 
such things as j-plane representation of the scattering amplitudqfj or its saddle point value 
do not have to be invoked in deriving the ratio (11571) . But, we still find the expression 
interesting, in that not only p and Acs. but also 7efr. and B are predicted here to depend 
on the kinematical variables (x, t,g^) through a single value j*. Furthermore, the limit of 
applicability of (I157p . the condition (I158p . is translated into 7cfr. < 0, which is to say that 



The real part to imaginary part ratio is often denoted by r]{t), but 77 is reserved for skewedness in this 
article. 

Dispersion relation 

Rc Ai-Hs,t) = fds''^^^^^^^ (159) 
TT y S' — S 

in the convolution form becomes a simple product form after Mcllin transformation, [Re (j, t)] — 
— cot(7rj72) X [Im t)]. This is why the ratio is better described in the j-plane language. 
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,2'iv'l + \A 



, 2j\/l + 



•2j 



-2i 



-27;/i + 



2jv'l + \/A 



(b) 



Figure 12: The integration contour in the z/-plane. The original contour (for t < 0) is on the 
real axis from negative infinity to positive infinity. The contour can be chosen as it is so long 
as iu* < 2, as shown in (a). When the saddle point is at iu* > 2, however, one will deform 
the contour to make it pass the saddle point u = u*; the process of deformation leaves an 
extra contribution coming from the pole of [1/ sm{7ij^)] at iu = 2 {jy = 2), as in (b). The 

jy = 4 pole corresponds to = ±2i^/ 1 + a/A in (b); when the saddle point value of iu* is 
even greater than that, the whole amplitude consists of contributions from the = 2,4 poles 
as well as from the contour passing through the saddle point [TT] . 



the GPD at that g^) increases in the DGLAP evolution. This observation is based 
simply on the j-plane representation of the scattering amplitude and an assumption that 
the kinematical variables {x,t,q^) are in the saddle point phase, and thus, does not rely on 
details of the hard wall model. 

The condition fllSSp is not satisfied, however, for sufficiently large g^; that is when '-/cs. > 
0, and GPD/PDF decreases under the DGLAP evolution. In this case, the contour of u 
integration should be taken as in Figure [T2] (b) , and the structure functions are expressed as 
a sum of contributions from j = 2 {iu = 2) pole and a continuous integration whose contour 
passes through the saddle point. For even larger g^, the saddle point value iu* may be as 
large as A^/^ and j* as large as 4, 6, ■ ■ ■ , and the structure functions are given by a sum of 
contributions from some finite number of poles j = 2, 4, 6, ■ ■ ■ < j* and the one from the 
saddle point approximation. Using the Kneser-Sommerfeld expansion of Bessel functions, 
these pole contributions (j = 2, 4, ■ ■ ■ < j*) can be written in the form 



ReZ 



m 



-7fehPn(r«j>)77*7*Pn("^in)- 



(160) 
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Therefore, they are interpreted as t-channel exchange of the n-th Kaluza-Klein mode of 
the spin j < j* state in the graviton trajectory [TT]. The remaining continuous integration 
should then be regardecj^^ as the w-channel exchange of j < j* states and all the sphere level 
contributions associated with j* < j states. 

Let us first focus on the < — t < region. Then, with a simple argument like we had 
in section |4.1.3|, one can see that the j = 2 {iv = 2) pole contribution dominates, other poles 
give rise to subleading corrections, and the saddle point contribution is even smaller than 
them [n]. 

The real parts of the structure functions are well approximated by the saddle point con- 
tribution for small enough g^, but they are expressed as a sum of j = 2,4, ■■ ■ < j* pole 
contributions and the saddle point contribution, when becomes large enough to satisfy 
j* > 2, that is, 

^Jn(gM)^ = > 2. (161) 
In 



This transition associated with the change in is understood as follows; it is usually better 
to treat the Pomeron trajectory exchange as a whole, not as a sum of exchange of individual 
spin j = 2, 4, ■ ■ ■ particles, because the amplitude of a spin j particle exchange has an ever- 
increasing factor Large virtuality of the photon and the photon wavefunction localized 
near the UV boundary, however, introduces 7*-7*-[spin j string] coupling suppressed by 
powers of (A/g), and the amplitudes of higher spin j exchange are suppressed significantly. 
The transition means that, for large g^, the power suppression (A/g)'^*^-'^ becomes more im- 
portant than {l/xY for higher spin j stringy states propagating in the sphere amplitude, and 
the contributions from smaller j, that is, j = 2,4, etc. "stand out" from all the rest. That 
is the physical meaning of the transition observed in the rearrangement of the u integration 
contour. Therefore, we refer to this phase as 

• KK-sum Low-spin Phase or Spin-2 Phase, where the saddle point value is larger 
than 2 (which also means 7efr. > 0). 

The discussion so far can be extended to the region A^ ^ (— t) ^ g^; the discussion in 
section 14.1.31 can be applicable almost in its form, except that the factor (A/q)"'^^^ in fl9T]) 
needs to be replaced by {^/—i/q)''^^^ for A^ <^ — t; the saddle point j* = j„* is the one using 
( 195|) in the hard wall model. The spin j pole contributions can also be evaluated explicitly; 



In this article, we used the dilaton-graviton scattering as the gravity dual model of the DDVCS process. 
If a D-brane is used as the model of the target hadron (baryon), however, the scattering amplitude is not a 
sphere amplitude with four NS-NS vertex operator insertion. 
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they are 



7{J) / A ^ 2A-2-7(i) 

T\j/2) \aVXxJ \q) 



For a given value of x and (so that j* > 2), the real part of the structure functions 
scale as because the j = 2 pole contribution dominates. For sufficiently large 

momentum transfer (so that j* < 2), however, the real part is better described by the saddle 
point contribution alone, and the real part shows the same scaling behavior in as the 
imaginary part of the structure functions. Therefore, the real part to imaginary part ratio is 
given by 



p{t,x,r] = 0,q^ 



f I j for - t< A2 if ([161]) holds true, 

(±) m A= « while VA,^^=Mg^ » 1. 

^ given by (11571) otherwise. 

(163) 

For a given value of and t, the ratio in the small x limit is eventually given by (I157p or 
possibly by the one with j* replaced by ap,n(^)- 

4.3.2 Impact Parameter Dependence 

Although the impact parameter space description of the real part does not have such an 
interpretation as transverse distribution of partons, it still is important in discussing how the 
unitarity limit is reached in 7* + /i scattering. Since the unitarity limit is not achieved in the 
Bjorken limit A^ <^ q^ without an extremely small x (and extremely large W'^ ~ q^ /x), the 
following discussion may be only of academic interest, but we present the result anyway. 
The phase shift x the impact parameter space (6, z, z') is given by [TU] 

2 2A{z) 2A{z') 

X'^'^'^^s, 6, z') ~ /C(., 6, z, z') (164) 

s 

at the leading order 'ml/Nc expansion (sphere-level); the unitarity limit of the photon-hadron 
scattering in a "partial wave" h is indicated when | t^'^p^'"''' | ~ (9(1) for z ^ 1/q and z' ~ 1/A 
[ini [m [I2j • The saddle point in v integration in the Pomeron kernel is determined regardless 
of whether we are studying impact-parameter dependent structure functions or phase shift; 
for 1 <^ (^A), where only the Pomeron trajectory j = ap^i{t) of the lowest KK 4D hadrons is 
relevant, we find 



sphere, , n, l_[l+e-^*]/ 1 V*" W A \ 2+^0*) _^(.*)^ 
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as long as the saddle point is not too large, iu* < 2 and j* < 2. Here, both the imagi- 
nary (absorptive) part and real (diffractive) part are included. The phase shift above shows 
expected dependence on A^^, x and g^, because 1 < j* and 2 + ■y{j*) > 2 — jo > 0. 

When the saddle point becomes large, iu* > 2, and hence j* > 2, the diffractive (real) 
part of the phase shift needs to be treated separately, as we have already seen in studying 
the momentum transfer dependence. In terms of kinematical variables, iu* > 2 for the lowest 
Kaluza-Klein mode (n = 1) is equivalent to 

1 / dj^^n=l 



ln([g/A]/[xv^])/v^ \ 



ln(g/A) + 



bA\ >2. 



(166) 



For either large ln(g/A) or large (6A), the diffractive part is written as a sum of j = 2 pole 
contribution and the saddle point contribution, possibly with some other finite number of 
j = 4, 6, ■ ■ ■ pole contributions. The j = 2 pole contribution is always the dominant one as 
long as iu* > 2, and behaves as 



z~l/q,z'~l/A 



1 



A2 



X e 



-•"12 lb 



(167) 



where m2,i is the mass eigenvalue of the lowest Kaluza-Klein mode of graviton in the warped 
throat. This impact parameter dependent profile (1167P is also a j* — )■ 2 (and iu* — )■ 2) limit 
of the real part of fll65p . All these results are shown schematically in the phase diagram in 
Figure [T3l 

The total cross section a'' ^ can be estimated by using the the sphere-level phase shift 
^sphere high energy, provided that the absolute value of phase shift \x\ is of order unity or 
larger within the critical radius be characterized by |x'^^'^'^'^'^(6c)| ~ C^(l)- Thus, the asymptotic 
form of the total cross section becomes 

A2 



7*h 
C^tot 



1 



(168) 



as long as iv*{hc) < 2. For even smaller x (larger W) with a given g^, however, h^. and iv*{hc) 
become larger, and the condition iu*{bc) < 2 will eventually be violated. The asymptotic 
behavior of cr^ot^i^, 9^) then turns into 



a. 



7*h 
tot 



In' 



m 



2,1 
hh' 



1 W'A 



2\2-\ 



(169) 



The total cross section of two light hadrons a (s) can also be discussed in a similar way; it 
can be calculated by replacing with A^. Then, the asymptotic behavior of a 



hh' I 



hh' 
'tot 



1 



In' 



(6c)' 



1 
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A2 



s) becomes 
(170) 
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Figure 13: A schematic picture describing qualitatively different behavior of Im x{b) (panel 
(a)) and Re x{b) (panel (b)) as functions of the impact parameter b. Behavior is different for 

different range of b as well as for different value of g^, as shown in this figure. (I): e""^! ''^ 
(see fll53p for a more precise expression), (II): e M{9/A)/VAa-) ^ (III):e 2in{{9/A)/x/Aa-) ^ ^^ci (IV): 



for iu*{bc) < 2. In the very large s region {iv*{bc) > 2), it becomes the form which is already 
shown in [T2] : 

''''2,1 

4.4 Structure Functions and GPD 

We have so far used holographic QCD to calculate the DDVCS amplitude and its gauge- 
invariant structure functions. While the scattering amplitude is all the necessary information 
in describing experimental data, a little more theoretical object GPD represents internal 
structure of hadrons more directly@ In this subsection, we will argue that GPD is a well- 
motivated notion separately from the scattering amplitude or structure functions even in 
strong coupling regime (in holographic QCD), and show how to extract GPD from holographic 
calculation. 

In the real world QCD, where only fermion partons (quarks and anti-quarks) are charged 

^°At least in weak coupling regime, GPD is also considered important, not just scattering amplitudes, also 
because GPD can be used to describe amplitudes of multiple different scattering processes. 



Ar2A2 



(171) 
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under the external probe (photon) of DDVCS, quark GPD in a scalar hadron h is defined by 
2 J ^^^''■(%2)|^(-2^)n-7V^(+2^)l%i)) = i^.(^,r/,t;g'), = 

(172) 

Its Mellin moment for an even integer j G 2N is 



poo 

/ dx x^"^ \iiq{x, 7], t; q^) + Hq{x, ?7, t; 
Jo 

/+00 i / ■ \ 

dxx^-'Hq{x,r],t;q^) = E^^^W ( ^ J {-2^)" = A,{r],t), (173) 

where Aj,fc(t) is the matrix element of the spin j (symmetrized traceless) twist-2 operator, 

fc=o ^''-•^ ^''^ o-e5j 

Here, ^ = (^ — ^)/2. The quark-anti-quark GPD, Hg{x,r],t) + Hg{x,ri,t), is regarded as 
the inverse Mellin transform of A^^\ri, t), a holomorphic function of j which becomes Aj{ri, t) 
for j e 2N. 

Since matrix elements of gauge-invariant local operators are well-defined notion even 
in gravitational dual descriptions, we characterize GPD in holographic calculation as the 
inverse Mellin transform of a holomorphic function of j, A^^\ri,t), which becomes matrix 
element of "twist-2" spin j operator for j G 2N. We have already seen in section 14.1.31 that 
[^hh¥*{j, t)]i/e (I120p is a holomorphic function of j which becomes Aj^k=o of some "twist-2" 
spin j operator for j G 2N. Thus, we can define a GPD renormalized at /i = 1/e as the 
inverse Mellin transform of [T(j,t)]^; (x,t) dependence of GPpRI can be calculated by using 
holographic QCD: 

H{x,v = 0.t;f,') = J ^^x-^[T{j,t)l. (175) 

The GPD obtained in this way satisfy the expected DGLAP evolution, as one can easily see 
from the e dependence of [T{j,t)]^^^^^ in fll2Up . 

Obviously GPD can be estimated in completely the same method as we have done in 
earlier sections for the structure functions. For the saddle point phase, 

i/(x,r/ = 0,t;/i2)~ (^-J l^-j [T{j\t)]A. (176) 



Discussion in this article needs to be refined in order to study skewedness dependence. 
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Note that the saddle point value now depends on /i instead of q. The form factor F for GPD 
defined in ( 199!) is the same as ( llOOp obtained by simply taking the ratio of structure functions 
at arbitrary t to the ones at t = 0. Therefore, all the discussion on the form factor at the 
end of section 14.1.21 can be read as statements on the form factor of GPD in the saddle point 
phase of strong coupling regime. Note that the (x, t, /i) dependent form factor fllOOD and the 
effective mass scale A^g are compatible with the renormalization group evolutiono For the 
leading pole phase, 

//(x., = o,,,^)~(-) (-) . (177) 

This is almost the same as the structure functions at the leading pole phase flll2p . 

Certainly one can calculate a "GPD" in holographic QCD, but it is not clear at this 
moment GPD of which parton it is. We adopted an asymptotically conformal (AdS^) geom- 
etry for holographic calculation, which means that the theory has all the particle contents of 
some conformal gauge theory. The reduced matrix element T{j,t) obtained in holographic 
calculation must be that of some linear combination of "twist-2" spin-j operators of gluons, 
fermions and scalars of the gauge theory. It must be possible, at least, to specify the linear 
combination that has the smallest anomalous dimension. Purely AdS^ background may also 
be replaced by a more realistic background. One could think of a couple of ways to refine 
the physical meaning of GPD obtained in this way, and possibly to obtain more results, but 
all of such improvements are beyond the scope of this article. 

Relation between the structure functions of the Compton tensor T^'^ and GPD in the 
strongly coupled regime (holographic calculation) remains quite similar to the one in the 
weakly coupled regime. To see this, let us first remind ourselves of relation between them in 
the weak coupling regime. Contributions from matrix elements of quark / gaugino twist-2 
operators to the Compton tensor T'"^ are formally written in the form of operator product 
expansion as 

^ Yl -4-^ E(c'..-);:r...p..r...JMP2)i(9™)--- [^7(^3)^-v]'"''^' impO); (i78) 

j=l m=0 

in perturbative QCD, the OPE coefficient (C'j,m)p^...pjcri---cr,„ can be calculated order by order 
in QCD coupling a^. Gauge-invariant structure functions Vi^...^^ can be extracted from the 



s^Such models of GPD as A^^^ = mlg in [671 and A^l ^ [a'(l - x)^ ln(l/x) + B(l - x)^ + Ax{l - x)] in 
[68j introduce ansatz of (x, t) dependent GPD profile at a given renormalization scale. The profile at other 
renormalization scale needs to be determined numerically by following the DGLAP evolution. 
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expression above. At tree-level in expansion, for example, the quark-ant i- quark contribu- 
tions to (some of) the structure functions at 77 = are given by 



i=2 



^ rAj,o{t), V^S^V ^ -- 

I ^ 2 x-> 



Aj-o(t)+4a;M,,2(t) 



(179) 



Thus, the structure functions at tree level at 77 = are expressed in terms of j-plane integrals 
(inverse Mellin transforms) as in 



Vi{x,7] = Q,t) 



e 2 



:i8o) 



l — x/^ + ie 1 + x/^ — ie 



gV3(x,r/ = 0,t) 



-1 



4i sin(7rj) x^ 

6 Jo ^2 2 
5\ 1 



J +2 



j(j-l)9(2r/)2 



6\ „fl 



+ 1- 



H^) (ei,r/ = 0,t) + 



:i8r 



Since the non-skewed structure functions in the complex j-plane representation are given by 
products of the signature factor — [1 -l- e~'^^^]/ sin(7rj), operator reduced matrix element Aj 
(or its ?7-derivative) and its OPE coefficient Cj, the structure functions become convolution 
of the inverse Mellin transforms of those factors. The factors [{^ — x + ie)~^ + {!^ + x — ie)~^] in 
f ll80p and f ll8ip originate from the inverse Mellin transform of the signature factor [20] here, 
but it also arises from the light-cone singularity of the fermion parton in direct calculation 
of the Compton tensor using Feynman diagrams [91] . 

In the holographic calculation of the structure functions of the DDVCS amplitude, we 
have seen that all the structure functions Vi^...^^ are given by /o,i, and /o,i are given by 
V integration (!78|) . The //-integration becomes j-plane integration, when the integration 
variable is changed through j = j^,. The integrand is a product of the signature factor 
— [1 + e~'"'^^]/ sin(7rj), reduced matrix element [r(j, t)]^ of some "twist-2" spin j operator and 
a remaining j-dependent factor, which is to be interpreted as the OPE coefficient and (q-p)^ oc 
x~^ . Therefore, the structure functions of DDVCS amplitude are given even in the strong 



coupling regime by a convolution of a "propagator" [(^- 



-x+ie] 



'{C,+x—ie) GPD (inverse 



Mellin transform of [r(j, t)]^) and inverse Mellin transform of OPE coefficient, just like in 
the weak coupling regime. The OPE coefficients in the strongly coupled regime, however. 
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are calculated through the vertex operator OPE on world sheet; the OPE coefficients may 
be improved order by order in 1/ \/~X expansion, rather than in expansion in perturbative 
QCD. 

5 Implication for the Real World QCD 

We have so far studied DDVCS and GPD by using holographic QCD. In this section, we 
discuss how much we can apply our results in gravity dual to the real world QCD. We should 
keep in mind that the hard wall model or other UV conformal holographic models are not 
equivalent to the real world QCD. Dual gauge theories of UV conformal gravity descriptions 
remain strongly coupled, even in UV scale, if the gravity descriptions are to be reliable. If 
one considers an asymptotic free holographic model, so that it matches on to the asymptotic 
free QCD of the real world, then such a gravity description will not be a useful framework 
in the UV region of the holographic radius. Although there is such a difference, and a big 
gap that is very difficult to ffil, it is also worthwhile to keep in mind that there are several 
features in GPD and structure functions that are shared by gravity dual descriptions and the 
real world QCD / perturbative QCD. 

Factorization theorem holds for (D)DVCS and leptoproduction of mesons in the real world 
QCD [921 193]. The structure functions are given by a convolution of hard kernel and GPD 
as in ( I180[ I18ip . The hard kernel is calculable in the perturbation of the gauge coupling, 
whereas the GPD is incalculable, and should be given as non-perturbative input. It would 
be nice if GPD profile can be determined by experimental data alone, but that is known 
not to be possible. Even in setting constraints on the proffie by using data, GPD has to be 
modeled and parameterized, because observable scattering amplitudes involve convolution of 
GPD, not just pure GPD, and furthermore, depend on GPD only in the limited region of 
kinematical variables. It is thus inevitable to construct physically motivated model of GPD. 

Modeling and parametrization of GPD has already grown into a large field, as one can 
see in review articles (e.g., [25| [26]). The simplest model imaginable assumes factorized form 
of GPD into PDF g(x) and a t-dependent factor [Ml [95] . 

H{x, ri = 0,t) = q{x)F{t), (182) 

and models that mimic small x Regge behavior were also considered [051 ETJ [HSl IHH] : 

i7(x,r/ = 0,t) =g(x)e*^(^) (183) 

with some x-dependent function g{x). In order to implement basic theoretical requirements 
on GPD, such as skewedness-polynomiality and proper DGLAP evolution, however, it is 
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more convenient to have a framework of parameterizing GPD where such requirements are 
implemented systematically. There are two such frameworks of systematic parametrization. 
The first one is double distribution |100[ IIOH 1102] along with D-term [103] ; transformation 
form this parametrization to GPD (Radon transformation) and its inverse is known |104[ 
I105[ 1106] . GPD models in this parametrization framework are found, for example, in |104[ 
1107] . The alternative framework of systematic parametrization is the coUinear factorization 
approach, which is also known as dual-parametrization. In this approach, models of GPD are 
constructed as an amplitude on the conformal moment space |108l I109[ I110[ IllH 1112] : GPD 
is given by a transformation |111] that becomes inverse Mellin transformation of t)]^ for 
zero-skewedness. 

Lessons from gravity dual calculation of DDVCS amplitude and GPD can be passed on 
to the understanding of GPD of the real world best in the form of j-plane representation, 
[r{j,t)]fj_, that is, in the coUinear factorization approach (dual parametrization). The trans- 
formation between the j-plane representation and Bjorken x representation is just a pure 
mathematical one, and does not rely on perturbation theory or even on existence of a cal- 
culable theoretical framework. The description of the (D)DVCS amplitude in j-plane is also 
based on OPE, which once again does not rely on perturbation theory [2]. While there is no 
way calculating hadron matrix element [r(j, t)]^ in perturbation theory, or in lattice gauge 
theory for complex valued j, various gravity dual models yield at least some results (if not 
truly faithful to the real- world QCD), as we have seen in earlier sections. 

The j-plane representation of the scattering amplitude and GPD also plays an essential 
role in characterizing the following three phases of the behavior of GPD and (D)DVCS 
amplitude: 

• j* < Re ap^i(t): leading pole phase (or leading singularity phase), 

• Re ap^i(t) < j* < 2: saddle point phase, where 7cff. < 0, 

• Re ap,i(t) < 2 < j*: spin-2 phase, where < 7es.. 

The three phases are also interpreted as Low-KK Spin-sum phase, KK-sum Spin-sum phase, 
and KK-sum Low-spin phase, respectively, from above to below. The distinction between 
the last two phases is absent in GPD or in imaginary part of the scattering amplitude (at 
leading order in l/N^ expansion), and we simply refer to the last two phases as saddle point 
phase. Earlier sections of this article described such a phase structure by relying on gravity 
dual calculation, but the essence is in the relative position of the leading singularity, j = 2 
pole and the saddle point in the j-plane representation, not so much on details of gravity 
dual descriptions. As we have studied in sections 14.1.11 and 14.1.31 the parton dynamics 
characterized by 7cfr. and Acfr. and the phase structure Figure [9] in {q'^,x) plane (fixed t slice) 
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in gravity dual are surprisingly similar (qualitatively) to the one we know in the real world 
QCD. This similarity is traced back to the fact that the saddle point is determined by ( 193!) in 
perturbative QCD as well as in holographic QCD, and that the anomalous dimensions 7(7) 
in both frameworks are also qualitatively similar [3]; especially the two properties 

^>0, '^«> (184) 

are shared for certain range of j. 

Most of kinematical regions explored in DIS experiments should be in the saddle point 
phase (including the spin- 2 phase in D1S)|^ as is clear from the (x, g^) dependence of 7cfr. 
and AefT. in DIS |113]Fn Anomalous dimension 7(j) in the weak coupling regimj^ should be 
used, because of the value of Aefr. (depending weakly on q"^) closer to than to 1 |113] . The 
GPD, then, is approximated by the expression f ll76p for some form factor [r(j, t)]\ evaluated 
at j = j*- This suggests (if the skewedness dependence that we did not study in this article 
does not screw things up) that various parameters characterizing the DVCS differential cross 
section^ 



it A* H A J • <^«^' 

namely, 5 = 4(j* — 1) and n = '~f{j*) + 2j* above, as well as the real part to imaginary 
part ratio p and the t-slope parameter B, are all functions of the saddle point value j*. 
Because of the way the saddle point value is determined, f l93|) . all the four parameters are 
more sensitive to ln(g/A) than to ln(l/x) or ln(l^/A) in small x, for the reason that we have 
already seen in (11331) . All of 5 = 4Acfr., n and p are increasing functions of j*, and 

hence they increasd^l for larger ln(g^/A^). They should decrease for larger ln(l/a;) or larger 
\n{W/A), but only weakly. Although it is not immediately clear whether the slope parameter 
B increases or decreases for larger j* (and hence for larger ln(g/A)), yet we found by using 
the spin j form factor in the gravity dual hard wall model (in section I4.1.4p that the slope 



^•^The leading singularity phase tends to be realized in smaller x, lower and less negative t, and might 
be found experimentally. It is difficult, though, to tell theoretically where the phase boundary is located in 
the [x, t, q^) variable space in a theory with running coupling constant. 

^""^It is not very clear in low q^ region, however, whether Acff. increases for larger [114) . c.f. Figure [5] (a). 
The leading pole phase may be hidden there. 

^^There are various approximation schemes in perturbative QCD, such as DLLA, BFKL and double leading 
approximation. 

^^This parametrization of DVCS cross section by ((5, n) (conventionally adopted in data analysis) corre- 
sponds to /, - (W^/A)2+V2(A2/g2)«/2^ 

67 The HERA measurements give 5 = 0.44 ± 0.19 for q^ = 2.4 GeV^, 6 = 0.52 ± 0.09 for q^ = 3.2 GeV^, 
,5 = 0.75 ±0.17 for q^ = 6.2 GeV^, 5 = 0.84 ±0.18 for q^ = 9.9 GeV^, and 5 = 0.76 ±0.22 for q^ = 18 GeV^ 
in ZEUS [ng, and 5 = 0.61 ± 0.10 ± 0.15 for q^ = 8 GeV^, ,5 = 0.61 ± 0.13 ± 0.15 for q^ = 15.5 GeV^, and 
,5 = 0.90 ± 0.36 ± 0.27 for q^ = 25 GeV^, in HI [78]. 



63 



parameter decreases for large g^, which is in nice agreement with the HERA measurements 

[761 Eg EH] H 

It is also worthwhile to remind ourselves that GPD models in the coUinear factorization 
approach (dual parametrization) in the saddle point phase are always compatible with the 
renormalization group evolution. In the approximate form (I176p . the renormalization scale 
H dependence comes into the form factor [T{j*,t)]\ through the saddle point value j*. 

All of these phenomenological as well as theoretical successes of the GPD models fll76p 
come directly from the nature of the saddle point phase, as we have already seen above. One 
more necessary ingredient is that the slope parameter 



2|ln[r(j,t)] 



= B (186) 

-t~[0~l]GcV2 



is a decreasing function of j; this property of the hard-wall model prediction was used in 
explaining the dependence of the t-slope parameter. It is a blessing indeed that the 
various features of experimental data are attributed to such a small number of conditions. 
To our knowledge, the first model of GPD in this category was that of |116] . which introduces 
an ansatz 

m t)]A - m = fi - ^ V" ; (187) 



j-a(t)' V 

a and n are parameters of the model. All possible holographic models will also have their 
own predictions of the reduced matrix element (spin j form factor) [r(j, t)], and hence suc- 
cessful models of GPD in this category are obtained from these models, as long as (11861) is a 
decreasing function of j. 

A general lesson from gravity dual calculations of the spin j form factor [r(j, t)] is that 
even the leading trajectory (containing graviton) in 5-dimensions (or in 10-dimensions) be- 
comes a Kaluza-Klein tower of infinite Pomeron trajectories for hadron scattering in 4- 
dimensions@ The form factor is given in the following form 

[r(j,t)]A = E7^' (188) 
^ t - m^,„ 

n J'"- 



^^In the leading pole phase with an anomalous dimension 7(j) in the weak coupling regime, the slope 
parameter of the (D)DVCS differential cross section increases for larger ln(l/x), and depends only weakly on 
ln(g/A); 7(j) cx as is used in (|135p . This prediction does not fit very well with the HERA data [751 [771 [75] . 
In this article, we did not work out gravity dual prediction for the t-slope parameter B in the spin-2 pole 
phase, where 7off . > 0. The contribution from the real part dominated by a tower of spin-2 glueball exchange 
should be important there. 

^^On top of this tower of infinitely many Pomeron trajectories, gravity dual descriptions predict yet another 
tower structure of trajectories, because stringy states form a tower of daughter trajectories already in the 5D 
(or lOD) description. 



64 



where the n-th term corresponds to the n-th Kaluza-Klein excited Pomeron trajectory; t = 

„ is the (mass)^-spin relation. After summing up all the contributions in a Kaluza-Klein 
tower (1188p . however, the form factor has a power-law dependence in (— t). This property is 
understood simply by using Green function in warped 5-dimensions; details of the background 
geometry of gravity dual models are irrelevant. This power-law behavior at large (— t) nicely 
matches on to the power law behavior expected theoretically from naive quark power counting 
[SZl ESI ES EDI EI], albeit only at qualitative level. 

It should be kept in mind, however, that gravity dual descriptions are not an ideal and 
flawless framework in calculating the form factor. Although a gravity dual model with a 
carefully chosen background may serve as a faithful dual to the QCD of the real world, it 
should work as a reliable framework of calculation only in the IR region of the holographic 
radius. The UV region of the geometry should correspond to the real world QCD in the 
weak coupling regime, and hence the curvature of the background should be too large for 
a reliable calculation. It will be nice if it is possible to make even a crude estimate of the 
error in gravity dual calculation of the form factor (e.g., the power-law behavior from UV 
conformal gravity dual models is dictated by the conformal dimension associated with the 
target hadron, whereas the power is determined by the number of valence partons in the 
(empirically successful) naive power counting rule), and to use perturbative QCD in some 
way or other in controlling or reducing the error. 

Such a dream may not be totally unrealistic, because at least one can compare singularities 
in the j-plane representation. Gravity dual models for asymptotic free running coupling (at 
least up to some high energy scale), e.g., [7S], can be used to calculate the isolated poles 
of [r(j, t)]. It is just to solve a Schrodinger equation [3]. Poles in the large Re j region 
are reliable, because the corresponding wavefunctions are localized in the IR region of the 
warped spacetime, but those with smaller Re j are not. On the other hand, the BFKL theory 
can calculate a spectrum of its kernel. The spectrum is discrete, when the running coupling 
effect is taken into account, and the discrete spectrum is mapped into poles in the complex 
j-plane [29] . The spectrum of the BFKL theory in large Re j region, however, corresponds 
to wavefunctions that are dominantly in the small k_i_ (gluon transverse momentum) region, 
where the coupling is large and perturbation in does not work well. Thus, the BFKL 
theory might be used in determining the j-plane singularities of [r(j, t)] in a region (smaller 
Re j) where the results from gravity dual cannot be reliable. 

In the calculation of the form factor, the spectrum of the Pomeron kernel (poles in the 
j-plane) is not the only necessary information. We also need corresponding wavefunctions 
of the spectrum, which are to be multiplied by a hadron impact factor, and integrated over 
the holographic radius z, if we are to use the language of gravity description. Based on the 
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(a) (b) 

Figure 14: One can expect straightforward generalization of various results in this article for 
the following processes, (a): DVCS/DDVCS or timelike Compton scattering. The final state 
hadron h' is allowed here to be different from the initial one h. For example, it can be an 
excited mode of the initial state hadron. (b): Vector meson production. 

similarity between the k±_ factorization formula of perturbative QCD and the expression (!62|) 
of hadron scattering amplitudes in gravity dual, correspondence between the k± coordinate 
of the BFKL theory and the holographic radius has been suggested |117l I118[ [3]. Thus, 
the integration (and integrand) over the z coordinate in the UV region may be replaced by 
that of the BFKL theory over the k± coordinate, in a crude attempt at improving the form 
factor [r(j, t)] calculated entirely in gravity dual models. All of such attempts, however, are 
beyond the scope of this article. 

6 Discussion 

In this article, we used holographic QCD to study double deeply virtual Compton scattering 
(DDVCS) amplitude and generalized parton distribution (GPD). Analysis presented in this 
article, however, is restricted to the scattering amplitude with vanishing skewedness; in order 
to exploit holographic methods to study non-perturbative aspects of deeply virtual Compton 
scattering (DVCS), analysis in this article needs to be extended to cover the case with non- 
zero skewedness. This subject will be covered in a future publication |87j . 

We have presented in this article an improved conceptual understanding of Pomeron 
couplings and Pomeron form factor, and predicted phase transition (crossover, to be more 
precise) in DDVCS amplitude. Such results, however, are not just for DDVCS and DVCS, 
but are more general in nature; some of the statements in this article hold true for some 
other high energy scattering processes only with minor (and almost obvious) modifications. 
An obvious generalization is to relax the constraint that the final state hadron h' is the same 
as the initial state hadron h] 7*/i — 'y^*^h' (Figure [T^ (a)). When the final state hadron h' 
and the initial state hadron h are different Kaluza-Klein modes of the same field in the 5D 
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effective theory, Pomerons (with vacuum quantum number) can couple to the transition from 
h to h'. One only needs to replace the impact factor Phh{z') by Phh'{z') = [^n{z')^n'{z')] 
in writing down the amplitude. Scattering amplitude of vector meson production 'j*h ^ Vh' 
(Figure UM (b)) and timelike Compton scattering 7*/i — 7*/i' (Figure UM (a) where the final 
state photon has timelike virtuality) can also be obtained by simply replacing the final state 
photon wavefunction properly. Since all of those processes involve an initial state spacelike 
photon with large virtuality, the initial state photon wavefunction localized strongly toward 
UV region gives rise to much the same qf and W"^ dependence of the Pomeron contribution to 
those scattering amplitudes. All of those processes involve a hadron-hadron-Pomeron vertex 
(lower blobs in Figure [H] (a, b)), and the non-perturbative vertex (form factor) in those 
processes should be essentially the same as in DDVCS. When the saddle point j = j*{x, g^, t) 
of the complex j-plane representation of the scattering amplitudes has a larger real part 
than all the singularities in the j-plane, the form factor depends on various kinematical 
variables (such as and g^) only through the saddle point value j*, and the form factor 
becomes that of "spin j* current". The t-slope parameter can also be calculated for these 
processes, but it will be obvious that the our argument in section 14.1.41 is so generic that its 
conclusion holds true for vector meson production and timelike Compton scattering as well. 
There is sort of universality in the Pomeron-hadron (/i)-hadron [h') form factor (coupling), 
regardless of whether it is used for DDVCS, DVCS, time-like Compton scattering or vector 
meson production (Figure [H] (a, b)). 

We have also seen in this article that the (Pomeron contribution to the) DDVCS ampli- 
tude shows three qualitatively different behaviors; the three phases are characterized by the 
position of the saddle point j* of the j-plane representation of the amplitude, relatively to 
the leading singularity and j = 2 pole of the signature factor, as summarized in section O 
The essence of this phase classification is shared also by DDVCS/DVCS with h' ^ h, vector 
meson production and timelike Compton scattering. Since all of these processes involve the 
kinematical variable qf, spacelike virtuality of the initial state photon, it is not surprising to 
see a crossover behavior in the differential cross sections of these processes. 
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